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ABSTRACT 

The  operation  of  a  two-phase  charge-coupled  shift  register 
is  discussed  and  analyzed.  A  parametric  study  for  simultane¬ 
ous  operation  with  more  than  one  input  signal  is  presented. 

The  output  circuit  is  analyzed  for  its  contribution  to  non¬ 
linearity.  An  analog  signal  processing  linear  model  of  the 
charge-coupled  device  is  proposed. 
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I .  INTRODUCTION 


Charge-coupled  semiconductor  devices  (CCD's)  were  inven¬ 
ted  in  1970  at  the  Bell  Laboratories  [Ref.  1].  Since  their 
invention  and  initial  experimental  investigation,  they  have 
held  forth  promise  of  low  noise,  high  dynamic  range,  analog 
delay  lines.  Great  effort  has  been  successfully  expended  to 
develop  charge-coupled  devices  into  self-scanned  solid-state 
image  sensors.  The  devices  are  held  to  be  attractive  for 
analog  and  digital  signal  processing  and  digital  memory.  The 
simpler  fabrication  of  charge-coupled  as  opposed  to  Metal 
Oxide  Semiconductor  Field  Effect  Transistor  (MOSFET)  or  bipo¬ 
lar  transistor  technologies  promises  increased  production 
yields  to  manufacturers.  The  net  result  is  a  quickening 
interest  in  charge-coupled  devices  for  signal  processing. 

A.  CHARGE- COUPLED  SEMICONDUCTOR  DEVICES 

Charge-coupled  semiconductor  devices  consist  of  closely 
spaced  MOS  capacitors  pulsed  into  deep  depletion  by  the  clock 
phase  voltages.  For  times  much  shorter  than  that  required 
to  form  an  inversion  layer  of  minority  carriers  by  thermal 
generation,  potential  wells  will  be  formed  at  the  silicon 
surface.  The  minority  carrier  charge  representing  the  infor¬ 
mation  will  be  stored  or  confined  in  these  potential  wells. 
The  propagation  of  the  information  is  accomplished  by  clock 
pulses  applied  to  the  electrodes  of  the  successive  MOS  capac- 
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itors  (i.e.,  charge-coupled  elements),  resulting  in  motion 
or  spilling  of  charges  from  the  potential  wells  that  become 
shallower  to  the  potential  wells  that  become  deeper.  Such 
propagation  of  signal  into  successive  minima  of  the  surface 
potential  produces  a  shift  register  for  analog  signals  having 
signal  transfer  efficiency  approaching  unity.  Such  analog 
shift  registers  may  be  used  for  various  signal  processing 
applications  such  as  electronically  variable  delays  or  self¬ 
scanning  photosensor  arrays. 

1 .  Multi-Phase  Charge-Coupled  Devices 

If  the  charge-coupled  structures  are  formed  with 
symmetrical  potential  wells,  at  least  three  phases  are 
required  to  determine  the  direction  of  the  signal  flow.  The 
use  of  more  than  three  phases  may  be  dictated  by  either  the 
construction  design  symmetry,  as  in  the  case  of  a  four-phase 
silicon  gate  overlapped  by  aluminum  structures,  or  special 
signal  coding  schemes  in  which  more  than  one  bit  may  be  pro¬ 
pagated  in  one  clock  cycle.  One  important  feature  of  the 
three-phase  system  is  that  it  may  be  used  for  a  bidirectional 
charge-coupled  channel  in  which  the  flow  of  information  may 
be  reversed  by  reversing  the  timing  of  two  of  the  three  phase 
clocks . 

2 .  Two-Phase  Operation 

Two-phase  operation  requires  that  the  charge-coupled 
structures  be  formed  so  that  the  potential  wells  induced  by 
the  phase  voltage  pulses  are  deeper  in  the  direction  of  the 
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signal  flow.  In  this  case,  as  one  phase  voltage  is  lowered, 
the  resulting  potential  barriers  force  a  unidirectional 
signal  flow.  It  should  be  noted  that  a  one-clock  operation 
can  also  be  obtained  in  a  two-phase  charge-coupled  shift 
register  if  a  proper  dc  voltage  is  applied  to  one  of  the 
phases  [Ref.  2].  Furthermore,  a  true  single-phase  or  a  uni- 
phase  charge-coupled  structure  [Ref.  3]  can  be  formed  by 
replacing  the  dc-biased  phase  with  a  structure  involving  a 
fixed  charge  in  the  oxide. 

B.  ANALOG  SIGNAL  PROCESSING  APPLICATIONS 

Greatly  overshadowed,  until  recently,  by  imager  develop¬ 
ment,  analog  signal  processing  is  a  fertile  area  for  charge- 
coupled  devices.  Electronically  variable  analog  delay  lines, 
correlators,  transversal  filters  and  recursive  filters  are 
among  the  possibilities.  It  is  assumed  in  all  of  these  appli 
cations  that  the  process  is  linear,  that  the  technique  of 
superposition  applies  to  the  input-output  relationships.  It 
is  this  assumption,  or  at  least  the  linearity  of  charge- 
coupled  devices,  that  is  examined  in  this  study. 

1 .  Transversal  Filters 

Filters  using  the  variable  delay  function  of  the 
charge-coupled  device  without  employing  signal  feedback  are 
classified  as  transversal  filters.  These  filters  are  charac¬ 
terized  by  a  finite  memory  and  little  or  no  phase  distortion. 
Sharp  cutoff  in  such  a  filter  can  only  be  achieved  with  a 
large  number  of  terms,  or  delays,  requiring  many  delay  units 
to  be  fabricated.  Such  filters  inherently  require  big  charge 
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coupled  device  arrays. 

2 .  Recursive  Filters 

Recursive  filters  employ  negative  feedback  of  the 
output  of  the  circuit  to  the  input.  Such  filters  enjoy 
infinite  memory  and  use  fewer  terms  than  do  transversal  fil¬ 
ters.  The  economy  of  terms,  and  hence  delay  lines,  usually 
results  in  poor  phase  characteristics.  It  is  simpler  to 
obtain  sharp  cutoff  using  recursive  techniques  than  with 
transversal.  This  represents  an  economy  in  charge-coupled 
structures  required. 

3 .  Variable  Delay  Lines 

Since  the  charge-coupled  device  moves  signals  along 
the  surface  of  a  channel  beneath  clocking  gates,  it  is  the 
frequency  of  the  clocking  signal  and  the  number  of  charge- 
coupled  elements  that  determine  the  delay  to  which  a  signal 
is  subjected.  Variation  of  the  clocking  frequency  can  be 
easily  achieved.  Such  variation  permits  simple  electronic 
control  of  the  delay  provided  by  the  device.  Use  of  low 
clocking  frequencies  provides  the  capability  for  small  long 
time  delay  units,  as  well  as  providing  a  fine  tuning  capa¬ 
bility. 
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II.  NATURE  OF  PROBLEM 


The  principle  of  superposition  is  so  pervasive  in  the 
area  of  analog  signal  processing  that  nearly  all  circuit 
design  and  analysis  is  dependent  upon  the  linear  nature  of 
the  components  used.  Specifically,  in  the  design  of  both 
transversal  and  recursive  filters,  the  process  is  assumed  to 
be  linear  as  a  fundamental  concept.  In  analog  signal  proc¬ 
essing,  it  is  vital  that  the  charge-coupled  device  be  oper¬ 
ated  in  such  a  manner  that  the  assumption  of  linearity  be 
valid.  How  to  assure  linear  operation  is  the  nature  of  the 
problem. 

A.  SAMPLED  NATURE  OF  THE  OUTPUT 

In  charge-coupled  devices,  an  output  signal  is  only  avail¬ 
able  when  the  clock  phases  have  moved  a  quantity  of  charge  to 
the  output  of  the  device.  This  quantity  of  charge  may  be 
proportional  to  an  input  signal  level  present  at  the  input  of 
the  device  at  a  time  in  the  past  equal  to  the  number  of 
charge-coupled  structures  divided  by  twice  the  clocking  freq¬ 
uency.  The  output  of  the  device  appears  as  a  sampled  version 
of  the  input  waveform  delayed  in  time  as  shown  in  Figure  1. 

The  sample  rate  and  the  delay  time  are  set  by  the  clock  freq¬ 
uency.  While  a  sampled  output  signal  can  hardly  be  consid¬ 
ered  linear;  provided  that  the  clock  frequency  is  high  rela¬ 
tive  to  the  signal  frequency,  the  linear  nature  may  be 
easily  provided  by  the  use  of  a  sample  and  hold  circuit,  or 
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in  some  cases  by  a  low  pass  filter. 


B.  ESTABLISHMENT  OF  THE  OPERATING  POINT 

The  experimental  devices  used  in  this  study  required  the 
setting  of  ten  variables  for  operation.  In  most  cases,  these 
variables  were  tuned  for  operation  only.  The  establishment 
of  an  operating  point  for  linear  operation  was  primarily  a 
function  of  the  dc-bias  levels  placed  on  the  input  ports  of 
the  device.  The  absolute  level  of  these  biases  was  effected 
by  some  of  the  other  parameters,  but  the  relative  value  of 
one  to  the  others  established  the  linear  region  desired.  An 
extensive  part  of  this  study  was  the  determination  of  the' 
bias  levels  necessary  for  linear  operation. 
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C.  TWO -PORT  MODEL 


When  electronic  devices  are  operated  so  that  non-linear 
distortion  is  small,  they  can  be  approximated  by  linear 
equations  and  circuit  elements.  In  such  cases,  all  the  pro¬ 
cedures  of  linear  circuit  analysis  may  be  applied  and  complex 
circuits  may  be  analyzed  with  a  minimum  of  effort.  Such  an 
approximation  is  a  linear  model  or  equivalent  circuit  of  the 
device.  This  model  can  be  used  to  determine  the  effect  of 
the  electronic  device  on  an  external  circuit.  It  is  import¬ 
ant  to  realize  that  a  model  does  not  represent  the  actual 
internal  behavior  of  the  device.  Most  active  devices  used 
in  electronic  circuits  may  be  modeled  as  a  two-port  linear 
system.  It  is  reasonable  to  assume  that  a  charge-coupled 
device  may  also  be  so  modeled. 

D .  OUTPUT  LINEARITY 

The  representation  of  the  signal  by  a  quantity  of  charge 
in  a  potential  well  within  the  semiconductor  produces  a  high 
impedance  sensing  problem.  The  output  circuitry  must  provide 
for  sensing  the  quantity  of  charge  on  each  sampling  by  the 
clock  phase,  provide  power  to  the  succeeding  circuits  with¬ 
out  diminishing  the  signal  charge,  and  finally,  it  must  dis¬ 
card  the  remnants  of  the  signal  charge  after  sampling.  If 
the  device  is  to  be  linear,  then  this  output  circuit  must 
also  be  linear. 


15 


III.  TWO-PHASE  CHARGE -COUPLED  SHIFT  REGISTER 


A.  TYPES  OF  TWO-PHASE  CCD'S 

The  charge-coupled  devices  initially  described  by 
Boyle  and  Smith  [Ref.  1]  require  three  or  more  phase  clocks 
to  obtain  directionality  of  signal  flow.  However,  for 
most  applications,  high  packing  density  and  better  per¬ 
formance  may  be  achieved  with  two-phase  charge-coupled 
structures.  As  illustrated  in  Figure  2,  the  asymmetrical 
potential  wells  or  barriers  in  the  surface  potential, 
needed  to  provide  directionality  of  the  information  flow 
for  the  two-phase  operation  can  be  achieved  by  incorpor¬ 
ating  into  the  charge-coupled  structures  one  of  the  fol- 
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lowing  features : 

1.  two  thicknesses  of  the  channel  oxide, 

2.  dc  offset  voltage  between  two  adjacent  gates  power¬ 
ed  by  the  same  phase  voltage, 

3.  two  levels  of  fixed  charge  in  the  channel  oxide,  or 

4.  ion-implanted  barriers. 

The  first  three  of  the  above  two-phase  charge-coupled 
shift  registers  can  be  conveniently  implemented  by  self-align 
ed,  closely  spaced  structures  in  the  form  of  polysilicon 
gates  overlapped  by  aluminum  gates  [Ref.  4] .  In  the  experi¬ 
mental  part  of  this  study,  a  charge-coupled  device  of  the 
first  type  above  was  used.  Described  in  this  section  in 
greater  detail,  it  was  a  two-phase  CCD  with  two  thicknesses 
of  channel  oxide  using  polysilicon  and  aluminum  gates. 

B.  DESCRIPTION  OF  EXPERIMENTAL  DEVICE 

The  device  used  in  this  study  was  fabricated  by  TRW, 
Incorporated,  for  the  purpose  of  experimental  studies.  The 
substrate  used  was  1.5  ohm- cm. ,  n-type  silicon  with  <111> 
orientation.  Construction  of  the  device  is  illustrated  in 
cross-section  with  a  plan  view  of  the  final  metallization  in 
Figure  3.  Six  p+-dif fusions  were  inserted  to  act  as  input 
and  output  diodes,  and  as  the  source  and  drain  diffusions  for 
two  output  circuit  MOSFETs.  Thermally  grown  silicon  dioxide, 
SiO^,  was  grown  to  approximately  1000  Angstrom  (A)  to  act  as 
channel  oxide.  Polysilicon  film  was  deposited  and  defined 
into  the  polysilicon  gates.  A  second  layer  of  silicon 
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dioxide  was  thermally  grown  for  the  aluminum  gates  to  a  total 
thickness  of  3000  A.  The  device  structure  was  completed  by 
opening  contacts  through  the  silicon  dioxide  to  the  p+-dif- 
fusions  and  the  polysilicon  gates.  A  final  metallization  of 
aluminum  was  deposited  and  defined  into  gates,  contact  lands 
and  leads.  The  entire  structure  was  approximately  33  mils 
long  by  14  mils  wide.  The  shift  register  input  gates  con¬ 
sisted  of  two  aluminum  gates  and  a  polysilicon  gate,  each 
0.5  mils  long  by  2.0  mils  wide,  and  a  polysilicon  gate,  1.5 
mils  long  by  2.0  mils  wide.  The  channel  gates  were  each  0.5 
mils  long  by  2.0  mils  wide,  resulting  in  a  channel  length  of 
2.0  mils  per  shift  register  stage.  The  output  gate  was  0.5 
mils  long  by  2.0  mils  wide.  The  two  MOSFETs  were  0.6  mils 
long  by  0.4  mils  wide. 

As  shown  in  Figure  3,  the  device  had  an  input  stage  con¬ 
sisting  of  the  source  or  input  diode,  D-l,  the  input  gates, 

G-l,  G-2,  G-3  and  G-4.  Separate  electrical  access  had  also 
been  provided  to  the  polysilicon  and  aluminum  electrodes  of 
each  phase;  i.e.,  4. - 1  (polysilicon),  $-1  (aluminum),  ^-2  (poly¬ 
silicon),  and  4-2  (aluminum) .  The  output  was  sensed  as  a 
current  flow  out  of  the  source  diffusion  or  output  diode, 

D-0.  The  output  gate,  G-0,  provided  a  means  of  setting  the 
surface  potential  in  the  vicinity  of  the  output  diode.  Two 
MOSFETs  were  diffused  onto  the  substrate,  with  one  having 
its  gate  and  drain  available  for  connection,  G-5  and  D-5. 

The  other  MOSFET  had  its  drain  and  source  available  as  D-6 
and  S-6.  The  second  MOSFET  gate,  the  source  of  the  other 
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MOSFET  and  the  output  diode  were  connected  together  on  the 
substrate,  but  were  not  available  off  the  substrate  for  ext¬ 
ernal  connection. 

The  device  was  mounted  in  a  forty  lead  flat  pack,  but  was 
not  sealed  or  protected  against  the  environment. 

C.  OPERATION  OF  THE  EXPERIMENTAL  DEVICE 

The  experimental  device  used  in  this  study  was  a  two- 
phase  charge-coupled  shift  register  employing  two  thicknesses 
of  channel  oxide  to  achieve  the  asymmetrical  potential  wells 
necessary  for  signal  flow. 

1 .  Charge-Transfer  Operations 

Assuming  essentially  zero  fixed  charge  in  the  channel 
oxide,  substrates  with  relatively  large  doping  concentrations 
are  required  to  obtain  a  substantial  difference  between  the 
surface  potential  under  the  polysilicon  gates  and  the  poten¬ 
tial  under  the  aluminum  gates  powered  by  the  same  phase  volt¬ 
age.  It  has  been  found  [Ref.  5]  that  the  potential  barrier 
formed  under  the  aluminum  gate  with  respect  to  the  surface 
potential  under  the  polysilicon  gate  of  the  same  phase  is  not 
constant  as  the  phase  voltage  changes .  The  maximum  amount  of 
charge  signal  that  can  be  stored  and  transferred  depends  in 
some  degree  on  the  frequency  and  the  waveshapes  of  the  phase 
voltages.  These  waveshapes  may  be  symetrical  with  equal 
rise  time  and  fall  time,  nonoverlapping,  or  overlapping.  In 
the  overlapping  case,  the  transfer  of  charge  is  preceded  by 
the  condition  in  which  both  phase  voltages  are  minimum.  If 
it  is  assumed  that  the  signal  charge  is  originally  contained 
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in  the  potential  well  beneath  the  phase-one  gate,  then  as 
the  phase  voltage,  is  changed  from  =  -23  volts  to 

V4, x  =  -13  volts,  the  surface  potential  under  the  phase-one 
gate  is  raised  to  a  higher  value,  the  charge  signal  is  trans¬ 
ferred  to  the  potential  well  under  the  phase-two  gate.  This 
mode  of  operation  for  two-phase  charge-coupled  devices  may 
be  referred  to  as  the  "complete  charge-transfer"  mode. 

The  experimental  device  used  in  this  study  was  operated 
in  the  complete  charge-transfer  mode  using  overlapping  wave¬ 
shapes  with  near  equal  rise  and  fall  times  for  the  two  phase 
voltages.  This  complete  charge-transfer  mode  of  operation 
is  illustrated  in  Figure  4  by  two  profiles  of  the  surface 
potential.  In  this  case,  at  time  t^,  the  charge  signal  is 
accumulating  in  the  potential  wells  under  the  phase-two 
gates;  and  in  the  second  half-cycle,  at  time  t2>  the  charge 
signal  will  be  transferred  to  the  potential  wells  under  the 
phase-one  gates. 

Another  mode  of  operation  of  two-phase  charge-coupled 
devices  can  be  achieved  if  the  surface  potential  under  the 
aluminum  gates  is  maintained  high  enough  in  all  cases  so  that 
the  potential  wells  under  the  polysilicon  gates  can  never  be 
completely  emptied.  In  other  words,  the  surface  potential 
under  the  polysilicon  gates  powered  by  the  minimum  phase 
voltage  is  lower  than  the  surface  potential  under  the  alumi¬ 
num  gates  powered  by  the  maximum  phase  voltage.  This  mode 
of  operation  may  be  considered  the  "bias  charge"  or  "bucket- 
brigade"  mode.  A  bias  charge,  or  background  charge  will 
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always  be  present  under  all  the  polysilicon  gates  in  this 
mode.  In  steady- state  operation,  the  bias  charge  will  be 
maintained  by  thermal  generation. 

The  transition  from  the  complete  charge-transfer  mode 
to  the  bias-charge  mode  of  operation  could  be  accomplished 
by  increasing  the  dc-bias  level  of  the  phase  voltages. 

The  presence  of  positive  charge  in  the  channel  oxide,  in 
the  case  of  n-channel  devices,  will  also  increase  the  bar¬ 
riers  under  the  aluminum  gates  sufficient  to  cause  bucket- 
brigade  mode  of  operation.  The  bias-charge  mode  of  operation 
of  two-phase  charge-coupled  shift  registers  is  very  similar 
to  the  operation  of  bucket-brigade  [Ref.  6]  shift  registers. 
In  the  bucket-brigade  device  case,  the  bias  charge  regions 
are  replaced  by  floating  diffusions;  otherwise,  the  oper¬ 
ation  of  the  two  types  of  shift  registers  is  very  similar. 

As  was  generally  known  in  connection  with  the  operation 
of  a  bucket-brigade  shift  register  [Ref.  6],  a  two-phase 
charge- coupled  shift  register  can  also  be  operated  with  one 
of  the  phases  dc-biased  at  a  level  used  on  the  other  phase. 
Since  no  switching  levels  are  used  on  the  dc  phase,  this  is 
referred  to  as  "one-phase"  operation.  This  one-phase  oper¬ 
ation  was  confirmed  with  the  experimental  device  used  in 
this  study. 

A  uniphase  charge-coupled  structure  requiring  only  one 
set  of  externally  controlled  gates  could  be  fabricated. 

Such  a  device  would  require  that  the  dc-biased  phase  gate  be 
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replaced  by  built-in  bias  in  the  structure.  This  bias  may  be 
obtained  by  the  intentional  presence  of  fixed  charge  in  the 
channel  oxide. 

2 .  Signal  Input  to  the  Experimental  Device 

Many  different  techniques  have  been  developed  for  the 
insertion  of  charge  in  charge-coupled  shift  registers.  For 
use  as  an  analog  signal  processing  device,  it  is  desirable 
that  the  technique  used  produce  a  charge  quantity  in  the  pot¬ 
ential  well  of  the  charge-coupled  shift  register  that  is  lin¬ 
early  related  to  the  signal  potential  or  current.  Further, 
the  technique  should  not  inject  noise  and  preferably  should 
demonstrate  high  transfer  efficiency. 

In  the  device  studied,  charge  injection  by  potential 
equilibrium  or  "sloshing"  technique  was  used.  This  method  of 
signal  charge  injection  is  demonstrated  in  Figure  5.  This 
diagram  shows  the  three  sequential  surface  potential  condi¬ 
tions  that  exist  during  one  transfer  cycle.  Also  shown  are 
the  relative  voltage  levels  required  to  produce  the  diagram¬ 
ed  surface  potentials  at  the  times  t^,  t2  and  t^.  The  input 
diode  is  pulsed  to  a  low  potential,  V^,  to  inject  charge 
into  the  potential  wells  formed  under  the  input  gates  Gl,  G2, 
G3  and  G4.  The  quantity  of  charge  injected  is  a  function  of 
the  level  to  which  the  diode  is  forward  biased,  the  duration 
of  the  pulse,  and  the  depth  of  the  potential  wells  under  the 
input  gates.  It  is  limited  by  the  magnitude  of  the  input 
pulse  and  the  potential  barrier  under  the  first  phase  gate. 
This  condition  is  demonstrated  by  the  surface  potential  dia- 
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At  the 


gram  in  Figure  5b,  the  condition  existing  at  time  t^. 
end  of  the  input  pulse,  the  input  diode  is  reverse  biased. 
Excess  charge  in  the  potential  wells  under  the  input  gates 
is  drained  through  the  reverse  biased  input  diode.  The 
charge  level  is  set  by  the  surface  potential  under  the  first 
input  gate,  Gl,  and  the  first  clocking  phase  gate.  The  quan¬ 
tity  of  charge  is  established  by  the  depths  of  the  wells 
formed  by  the  surface  potentials  under  the  input  gates,  G2, 

G3  and  G4.  This  condition  is  shown  in  Figure  5c,  existing  at 
time  t2’  Notice  that,  provided  that  the  input  pulse  produced 
sufficient  charge  to  fill  the  wells  under  the  input  gates, 
the  signal  charge  is  no  longer  a  function  of  the  input  pulse. 
A  major  source  of  input  non-linearity,  the  input  diode,  has 
been  thus  dispatched.  Following  the  disposal  of  excess 
charge  through  the  input  diode,  the  first  clocking  gates  are 
pulsed  to  a  low  potential  by  the  clock  voltage.  Accumulated 
charge  under  the  input  gates  is  now  transferred  to  the  poten¬ 
tial  well  xander  the  first  clocking  gates.  This  condition  is 
diagramed  in  Figure  5d,  as  the  condition  existing  at  time  t^. 
As  can  be  seen,  the  input  signal  charge  quantity  is  a  func¬ 
tion  of  the  surface  potentials  beneath  the  input  gates.  The 
determination  of  the  signal  levels  for  which  a  linear  input 
can  be  obtained  was  a  very  large  portion  of  this  study. 

3 .  Output  Circuit  Used  in  the  Study 

The  output  circuit  used  in  the  experimental  study 
employed  an  output  diode,  a  reset  MOSFET,  and  an  output  MOS- 
FET  connected  as  a  source  follower.  The  circuit  schematic  is 
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shown  in  Figure  6.  The  output  diode,  DO, is  reverse  biased  by 
the  drain  bias  voltage,  E^,  through  the  reset  transistor, 

Q5.  The  degree  of  reverse  bias  is  determined  by  the  surface 
potential  adjacent  to  the  diode  diffusion  in  the  charge-cou¬ 
pled  device  channel.  This  surface  potential  is  set  by  a  dc- 
bias  on  an  output  gate  not  shown  in  the  diagram.  This  output 
gate  is  an  aluminum  gate  immediately  adjacent  to  the  output 
diode  diffusion.  This  output  gate  could  be  connected  to  the 
phase  clock  voltage  of  the  gates  next  to  it,  but  experimental¬ 
ly,  better  noise  immunity  was  obtained  with  a  dc-bias.  When 
the  phase  clocked  gates  at  the  end  of  the  charge -coup led 
structure  are  sent  to  their  low  potential,  the  surface  poten¬ 
tial  well  developed  beneath  them  accepts  the  signal  charge  of 
minority  carriers.  The  presence  of  this  signal  charge  modi¬ 
fies  the  potential  at  the  cathode  of  the  output  diode.  Since 
the  charge  in  the  depletion  layer  of  the  diode  is  bound,  the 
potential  at  the  anode  rises  to  reflect  the  signal  charge. 

The  increasing  potential  at  the  anode  of  the  output  diode  is 
applied  to  the  gate  of  the  output  MOSFET,  Q6.  The  more  posi¬ 
tive  potential  at  the  gate  of  the  output  MOSFET  decreases  the 
drain  to  source  current  flowing  through  this  transistor.  The 
lessened  drain  current  results  in  a  more  positive  voltage  at 
the  output. 

At  the  completion  of  the  clock  phase  cycle  which  opened 
the  clocking  gates  at  the  end  of  the  shift  register,  the  sur¬ 
face  potential  barrier  again  blocks  the  transfer  of  charge 
signal  to  the  output  diode.  A  refreshing  pulse  is  applied 
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to  the  heretofore  cutoff  reset  MOSFET,  Q5 .  This  negative 
going  refresh  pulse  turns  Q5  on  and  restores  the  output  diode 
to  a  fully  reverse  biased  condition;  sets  the  quiescent  volt¬ 
age  on  the  gate  of  the  output  MOSFET,  Q6 ,  at  or  near  the 
value  of  the  drain  bias  supply;  and  provides  a  low  impedance 
path  for  the  previous  signal  charge  to  flow  to  ground  through 
the  reverse  biased  output  diode. 

This  output  circuit  permitted  isolation  of  the  small 
charge  packets  that  represented  samples  of  the  input  signal 
from  the  loading  effects  of  the  electronic  test  instrumenta¬ 
tion  that  was  used  in  the  study.  The  linearity,  that  is  the 
relationship  between  the  potential  produced  by  the  charge 
signal  at  the  cathode  of  the  output  diode  and  the  voltage 
output  of  the  source  follower,  of  this  circuit  will  be  ana¬ 
lyzed  further  in  another  section  of  this  study. 


Refresh  Pulse 


CCD  Channel 


D*0 


Oi 


OUTPUT  CIRCUIT 
Figure  6 
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IV.  EXPERIMENTAL  PROCEDURES 


The  experimental  study  of  the  two-phase  charge-coupled 
shift  register  reported  in  this  paper  evolved  into  six  areas: 

1.  To  accomplish  charge  transfer, 

2.  To  establish  linear  signal  transfer, 

3.  To  determine  the  device  gain, 

4.  To  determine  circuit  model  parameters, 

5.  To  establish  linear  operating  regions,  and 

6.  To  determine  the  characteristcs  of  the  output 
MOSFET. 

A.  TO  ACCOMPLISH  CHARGE  TRANSFER 

The  experimental  device  was  connected  as  shown  in  Figure 

7.  Using  parameter  values  provided  by  the  manufacturer  based 
on  a  computer  study  of  a  typical  device,  charge  transfer  was 
established  by  trial  and  error.  As  a  general  technique,  dc- 
bias  levels  were  applied  to  the  input  gates,  Gl,  G2,  G3  and 
G4.  Since  no  data  was  available  on  the  surface  potential- 
gate  voltage  relationship  of  these  gates,  these  levels  were 
determined  empirically  and  usually  were  found  to  be  the  most 
critical  of  all  parameters  effecting  the  performance  of  the 
device.  The  potential  on  the  output  gate,  Eq,  was  set  at  a 
previously  successful  value  of  -30  volts.  The  MOSFET  drain- 
bias,  E^Ej,  was  likewise  set  at  a  nominal  value  of  -25  volts. 
The  clocking  phase  voltages  were  set  at  the  manufacturer's 
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recommended  values.  The  dc-bias  level  and  pulse  amplitude  of 
the  refresh  pulse,  Vre£,  were  set  to  a  value  that  would  turn 
the  reset  MOSFET  off  and  on  (pulsing  from  +2  to  -1  volts  rel¬ 
ative  to  the  drain  bias) .  The  dc  level  and  pulse  amplitude 
of  the  input  diode  supply,  VDI,  was  then  adjusted  until  sig¬ 
nals  could  be  observed  by  oscilloscope  at  the  output  terminal. 
Minor  perturbations  were  made  on  all  parameters  to  maximize 
the  output  current  pulses.  A  typical  set  of  operating  wave¬ 
forms  are  shown  in  Figure  8.  For  these  waveforms,  the  dc-bias 
levels  were: 


1. 

MOSFET  Drain  Bias , 

EDD  " 

-25.0 

volts 

2. 

Output  Gate  Bias, 

E0  = 

-30.0 

volts 

3. 

Input  Gate  One  Bias , 

EG1  = 

-19.6 

volts 

4. 

Input  Gate  Two  Bias , 

eG2  = 

-20.0 

volts 

5. 

Input  Gate  Three  Bias , 

EG3  = 

-22.0 

volts 

6. 

Input  Gate  Four  Bias , 

EG4  = 

v«n 

B.  TO  ESTABLISH  LINEAR  SIGNAL  TRANSFER 

With  the  charge- coup led  device  providing  signal  charge 
transfer  as  above  without  an  external  signal  applied,  the  in¬ 
put  diode  pulse  amplitude  was  increased  slightly.  If  the 
output  pulses  showed  no  increase  in  amplitude,  then  potent¬ 
ial  equilibrium  input  had  been  obtained. .. the  signal  level 
was  not  a  function  of  the  input  diode  pulse.  Should  this  not 
occur,  the  input  gate  biases  were  modified  until  potential 
equilibrium  input  was  obtained.  A  small  ac-signal  level  was 
then  impressed  upon  any  of  the  input  gates  for  testing.  The 
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TYPICAL  PROCESSING  WAVEFORMS 

Figure  8 
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device  parameters  were  then  again  tuned  to  optimize  the  sig¬ 
nal  output  for  amplitude  and  linearity.  With  a  small  signal, 
linearity  was  usually  optimized  by  slightly  varying  the  dc- 
bias  level  of  the  input  signal  so  that  the  output  quiescent 
value  of  voltage  lay  halfway  between  the  no  signal  level  and 
the  fully  saturated  output.  The  phase  shift  between  the  in¬ 
put  and  the  output,  i.e.,  the  time  delay,  was  then  checked  to 
insure  that  shift  register  action  had  been  obtained.  The 
observed  delay  was  to  be  eight  times  the  period  of  the  clock¬ 
ing  voltage. 

C.  TO  DETERMINE  DEVICE  GAIN 

Once  the  device  was  set  into  an  operating  mode,  two  tech¬ 
niques  were  established  for  gain  measurement.  The  first 
termed  dc-gain  measurement,  involved  setting  the  dc-bias 
values  on  all  the  input  gates  save  one  at  independent  para¬ 
metric  values.  The  dc-bias  on  the  remaining,  or  now,  signal, 
input  gate  was  varied  at  small  voltage  intervals.  The  result 
ing  source  current  of  the  output  MOSFET  was  determined  and 
plotted  as  a  function  of  the  input  gate  potential.  The  sec¬ 
ond  technique  involved  inserting  a  small,  known,  ac  signal  on 
the  input  gate  and  measuring  the  output  signal  either  with  an 
oscilloscope  or  a  wave  analyzer.  The  small  signal  ac-gain 
thus  determined  by  the  ratio  of  these  ac  values  or  from  the 
slope  of  the  plot  of  the  first  technique  was  recorded  as  the 
gain.  In  the  actual  case,  the  first  technique  was  used  to 
the  exclusion  of  the  second  except  for  spot  checks  of  results 


33 


Use  of  the  peak  pulse  output  current  as  a  function  of  the 
input  signal  voltage  permitted  a  more  precise  determination 
of  the  device  gain.  The  slope  of  this  plot  multiplied  by  the 
value  of  the  output  resistance  yielded  the  small  signal  volt¬ 
age  gain. 

A  minicomputer,  an  ins trument/ computer  interface  unit  and 
two  digital  multimeters  were  used  for  semiautomatic  data  col¬ 
lection.  The  equipment  employed  was: 

Minicomputer,  TEKTRONIX  Model  31 
Calculator- Instrument  Interface,  TEKTRONIX  153 
Digital  Multimeter,  TEKTRONIX  DM  501 
One  voltmeter  was  connected  across  the  device  output  and 
the  other  was  connected  to  the  selected  input  gate.  A  series 
of  calibration  runs  were  conducted  to  determine  the  peak  out¬ 
put  pulse  current  as  a  function  of  the  output  voltage  read  on 
the  digital  voltmeter.  The  slope  of  this  curve  provided  a 
proportionality  constant  for  converting  the  dc  voltage  read 
on  the  digital  meter  to  a  peak  pulse  current  value.  An 
extremely  low  frequency  function  generator  was  connected  to 
the  selected  input  gate  as  the  signal  source.  A  ramp  func¬ 
tion  at  a  frequency  of  about  0.002  Hz  was  used  to  vary  the 
input  gate  signal  at  such  a  rate  that  the  input  voltage  would 
change  by  only  the  least  significant  digit  of  interest  in  one 
sample  period  of  the  digital  voltmeter.  The  computer,  which 
sampled  both  the  input  and  output  voltmeters,  was  programmed 
to  store  both  voltages  read  at  the  instant  that  the  input 
voltage  was  equal  to  a  comparison  value  stored  in  the  comput- 
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er.  When  the  read  condition  had  been  satisfied  and  the  data 
stored,  the  comparison  value  would  be  incremented  by  the 
delta  value  inserted  previously  by  the  operator.  The  com¬ 
puter  would  then  again  wait  for  a  matching  condition.  In 
such  a  manner,  it  was  posible  to  automatically  record  the 
input  and  output  voltages  for  as  many  values  as  the  operator 
designated  at  the  start  of  a  run,  commencing  at  a  given  input 
voltage  and  incremented  by  a  given  change  of  input.  After 
the  desired  number  of  data  points  had  been  recorded,  the  com¬ 
puter  would  print  out  the  values  of  the  input  voltage  and 
the  product  of  the  output  voltage  and  the  previously  deter¬ 
mined  proportionality  constant,  i.e.,  the  output  peak  pulse 
current . 

As  inputs  to  the  semiautomatic  routine,  the  operator 
need  only  type  in: 

1.  An  arbitrary  run  number, 

2.  Initial  input  voltage, 

3.  Desired  incremental  increase  in  input  voltage, 

4.  The  number  of  input  data  points  desired,  and 

5.  The  proportionality  constant  to  convert  output 
dc  voltage  to  peak  pulse  source  current. 

The  operator  would  then  start  the  function  generator  at  some 
value  of  output  less  than  the  inital  input  voltage  selected 
at  the  frequency  necessary  to  permit  the  voltage  to  change  at 
a  slow  enough  rate,  and  then  pursue  other  investigations 
until  the  computer  printed  out  the  data. 


35 


D.  TO  DETERMINE  CIRCUIT  MODEL  PARAMETERS 

Initially,  in  the  study,  it  was  felt  that  the  charge- cou¬ 
pled  shift  register  could  be  modeled  as  a  two-port  linear 
device.  To  support  such  a  model,  considerable  effort  was 
expended  in  efforts  to  determine  the  input  admittance,  out¬ 
put  impedance  and  transfer  functions.  In  fact,  none  of  the 
required  parameters  were  sucess fully  measured  except  for  the 
gain  as  described  above.  The  pulsed  nature  of  the  device, 
that  is,  the  presence  of  a  pulsed  waveform  at  the  output, 
even  with  no  input  signal  present  interferred  with  the  deter¬ 
mination  of  the  output  impedance.  The  many  signals  necessary 
to  cause  a  CCD  to  function,  even  at  a  quiescent  level,  pro¬ 
duced  a  need  for  many  connections  to  the  device.  These  con¬ 
nections  terminated  in  several  signal  generators  and  power 
supplies.  All  efforts  to  measure  the  input  admittance  of 
the  device  were  swamped  by  the  background  capacitance  and 
noise  of  the  experimental  setup. 

E.  TO  ESTABLISH  LINEAR  OPERATING  REGIONS 

With  the  multiple  input  gates,  the  possibility  of  using 
the  charge-coupled  shift  register  as  an  adder  or  subtractor 
of  analog  signals  is  obvious.  A  major  requirement  for  such 
operations  to  be  successful  would  be  that  the  device  be  lin¬ 
ear  with  respect  to  all  signals  inputed.  It  is  thus  neces¬ 
sary  to  know  for  what  values  of  input  voltages  on  all  gates 
will  linear  operation  result.  The  individual  dc  gain  curves 
were  evaluated  to  determine  the  regions  over  which  operation 
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was  linear.  These  curves  showed  a  prompt  departure  from 
linearity  at  a  high  and  a  low  value  for  whichever  input  gate 
was  being  used  as  the  controlling  gate.  These  departures 
occurred  when  the  potential  wells  were  nearly  filled,  satur¬ 
ation;  and  when  the  wells  were  nearly  empty,  that  is  the 
signal  charges  could  not  be  sensed  at  the  output. 


F.  TO  DETERMINE  THE  CHARACTERISTICS  OF  THE  OUTPUT  MOSFET 

At  the  request  of  the  experimenter,  the  device  manufact¬ 
urer  established  the  test  connections  shown  in  Figure  9  on 
the  output  MOSFET  of  four  devices  similar  to  those  used  in 
this  study.  Under  such  an  experimental  connection,  the  volt¬ 
age  between  the  gate  and  the  drain,  V^,  is  zero.  With  this 
potential  difference  at  zero,  the  gate  to  source  voltage,  Vq£ 
is  equal  to  the  drain  to  source  voltage,  V^g.  This  is  the 
simplest  arrangement  for  testing  a  MOSFET,  for  the  transis¬ 
tor  is  then  in  the  saturation  region.  From  the  basic  circuit 
theory  developed  about  the  Sah  MOSFET  model  [Ref.  7] ,  under 
conditions  of  saturation;  that  is,  when: 

Ivgs-vt^vdsI 


then : 


IDS 


f  \ 

y  e  e 

pro 

2t 

u 

OX  ' 

(VGS-V 


(1) 


where 


]x  =  average  surface  mobility  of  holes  in  the  channel 

t  =  thickness  of  the  oxide  over  the  channel 
ox 

er  =  relative  permittivity  of  the  oxide 

eQ  =  permittivity  of  free  space,  8.854x10  ^  farad/cm 
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SL  -  length  of  the  channel  in  the  direction  of  current 
flow 

W  =  width  of  the  channel 
Vgg=  gate  to  source  voltage 

VDS=  ^raan  to  source  voltage 
VT  =  threshold  voltage 

Collecting  some  of  the  more  or  less  constant  terms  from  the 
above  expression  (1) ,  permits  the  definition  of  two  new 
constants : 


1Jper£o 
k'  =  ; 

2tox 

(2) 

?|C 

II 

(3) 

l 

Use  of  these  new  consolidated  constants  permits  a  simplifi¬ 
cation  of  the  expression  for  the  source-drain  current  to: 


IDS  ’  -k<VGS-VT>2  (4) 

Note  that  the  parameter,  k,  is  a  design  variable  which  is 

a  function  of  the  ratio  of  the  channel  width,  W,  and  of  the 
length  of  the  channel  in  the  direction  of  current  flow,  l. 

The  parameter  k'  is  fixed  by  the  fabrication  process.  For 
purposes  of  nomenclature,  k  is  described  as  the  gain  para¬ 
meter,  and  k'  is  termed  the  conduction  factor. 

Using  the  same  transistor  model,  the  drain  to  source 
current  in  the  non-saturation  condition  becomes : 

IDS  "  -k[2(vGS-vT)vDS-(VDS)2]  (5) 


when: 


Ivgs-vtI 


>|v. 


DS1 


In  the  saturation  condition  depicted  in  Figure  9 ,  if 
drain  to  source  resistance  is  neglected,  a  plot  of  |  /lDS | 
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as  a  function  of  |vGS|,  (or  VDg  since  the  two  potentials  are 
equal  in  this  configuration)  yields  a  straight  line.  The 
slope  of  this  curve  is  /k ,  with  an  abscissa  intercept  of  |vT|. 
That  this  is  so,  may  be  quickly  established  by  taking  the 
square  root  of  both  sides  of  the  equation  expressing  the 
drain-source  current  relationship  to  the  gate-source  voltage 
in  saturation  (1) : 

fhs  ’  ■'r'<VGS-VT>  <6) 

Such  a  relationship  is  demonstrated  in  Figure  10. 

The  experimental  device  manufacturer  provided  the  source 
to  drain  current  and  source  to  drain  voltage  data  for  ten 
data  points  on  four  samples  of  the  output  MOSFET. 


MOSFET  TEST  SETUP 
Figure  9 
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V.  PRESENTATION  OF  DATA 


In  reviewing  the  data  obtained  in  this  study,  it  should 
be  remembered  that  it  represents  determinations  made  on  three 
devices  of  the  same  design;  and  as  near  as  possible,  the  same 
fabrication  techniques.  During  the  course  of  the  study,  no 
effort  was  made  to  produce  a  modified  device  that  would  be 
optimized  to  the  conditions  of  the  experiment.  Rather,  a 
stock  device  for  basic  study  of  charge- coupled  shift  register 
principles  was  used  to  study  analog  signal  processing. 


A.  DC  GAIN  MEASUREMENTS 


DC  gain  determinations  and  output  current  pulse  amplitude 
as  a  function  of  control  gate  voltage  measurements  were  made 
for  thirty- six  variations  of  the  independent  parameters.  In 
these  measurements ,  control  gate  4  was  connected  to  the 
clocking  signal  for  phase  one,  thus  becoming  part  of  the 
clocking  structure.  The  voltage  on  control  gate  one  was 
varied  as  the  independent  variable.  The  peak  pulse  current 
amplitude  of  the  output  signal  pulse  was  defined  as  the 
dependent  variable.  The  dc  potentials,  or  biases,  on  control 
gates  2  and  3  were  set  as  independent  parameters.  The  func¬ 
tional  relationship  demonstrated  in  Figure  11  is  typical  of 
the  data  obtained.  From  such  data,  the  device  gain  can  be 
obtained  by  multiplying  the  slope  of  the  curve  times  the  out¬ 
put  load  resistance,  as: 


„  _  dlc 
g  =  — RL 

dV 


(7) 
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OUTPUT  CURRENT  PULSE 
AMPLITUDE,  l0  ,  (pomps) 


-19.40  -19.50  -19.60  -19.70 

DC  VOLTAGE  ON  GATE  NO.  I,  Eq  (  ,  (Volts) 

Typical  input  characteristics ,  l0  -  EGJ 

Figure  II 
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or  g  =  96.74xl0_6x50xl03 
g  -  4.8 

A  gain  of  4.8  was  observed  over  most  of  the  useful  range  of 
input  signal.  At  the  limits  of  operability,  the  gain  was 
reduced. 

B.  LINEAR  REGIONS 

In  the  DC  gain  studies,  there  were  two  independent  para¬ 
meters,  the  voltages  on  control  gates  2  and  3.  Since  adjust¬ 
ment  of  these  parameters  changed  the  voltage  range  of  control 
gate  one  over  which  the  output  was  linearly  related  to  the 
input,  it  seemed  that  a  second  gate  could  be  used  to  permit  a 
dual  input.  In  such  an  arrangement,  input  signals  could  be 
applied  to  gates  one  and  two,  or  gates  one  and  three.  No 
study  was  made  to  determine  the  relationship  of  inputs  on 
gates  two  and  three,  or  for  inputs  on  three  or  more  gates  at 
the  same  time.  The  capability  to  respond  to  inputs  on  more 
than  two  gates  certainly  exists,  however. 

The  DC  gain  curves  were  studied  to  determine  at  what  volt¬ 
ages  on  gate  one  departure  from  linearity  would  occur  as  a 
function  of  the  voltages  on  gates  two  and  three.  This  para¬ 
metric  study  provided  the  regions  over  which  linear  operation 
could  be  acheived  as  a  function  of  the  voltage  on  two  control 
gates  with  the  voltage  on  the  third  gate  as  an  independent 
parameter.  The  resulting  loci  of  linear  operation  are  pre¬ 
sented  in  Figures  12  and  13.  The  range  of  voltage  used  on 
control  gate  one  was  kept  the  same  for  the  study  of  the  effect 
of  control  gate  two  and  the  effect  of  control  gate  three.  As 
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DC  VOLTAGE  ON  GATE  NO. 


Locus  of  linear  operation,  gates  I  and  2. 

F  i  gure  12 
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DC  VOLTAGE  ON  GATE  NO.  3,  E63,(  Volts) 


-17.00  -18.00  -19.00  -20.00 

DC  VOLTAGE  ON  GATE  NO.  I,  EG  j  ,  (Volts) 

Locus  of  linear  operation,  gates  I  and  3. 


Figure  13 
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can  be  seen  from  the  two  loci,  simultaneous  linear  operation 
can  be  achieved  for  signals  on  two  gates.  It  should  be  noted, 
however,  that  the  dynamic  range  of  simultaneous  inputs  is 
small,  approximately  0.1  to  0.2  volts. 

The  upper  limits  of  the  loci  were  not  defined  in  this 
study.  This  was  because  linear  operation  could  be  continued 
with  proper  adjustment  of  the  various  potentials  on  the  CCD 
as  the  signals  became  more  and  more  negative.  The  limit  of 
such  accommodation  would  appear  to  be  the  breakdown  potential 
of  the  various  gates  concerned.  A  "saturation"  effect  of  the 
potential  on  gate  one  can  be  seen  in  the  upper  limits  of  Fig¬ 
ure  13,  however. 

C.  AC  GAIN  CONFIRMATION 

Small  signal  AC  gain  determinations  were  made  on  the 
device  to  confirm  the  DC  measurements.  A  20  millivolt  (rms) 
signal  was  applied  to  control  gate  one  with  all  other  input 
gates  at  bias  levels  equivalent  to  those  used  in  the  DC  gain 
measurements.  The  output  voltage  was  measured  with  a  wave 
analyzer,  and  the  ratio  of  the  output  rms  voltage  at  the  sig¬ 
nal  fundamental  frequency  to  the  input  rms  voltage  was  record¬ 
ed  as  the  voltage  gain.  The  dc-bias  on  control  gate  one  was 
varied  as  an  independent  parameter  and  the  AC  gain  determined 
as  a  function  of  this  bias.  Figure  14  is  representative  of 
the  results  obtained.  No  significant  departure  from  this 
curve  was  noted  by  varying  the  input  signal  frequency.  The 
gain  level  of  4.75  observed  agrees  well  with  the  4.8  from  the 
DC  gain  determinations.  A  horizontal  translation  of  the 
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-19.4  -19.5  -19.6  -19.7 

DC  VOLTAGE  ON  GATE  NO.  I,  EG,  ,  (Volts) 

Figure  14 
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of  the  curve  along  the  abscissa  was  noted  when  the  independ¬ 
ent  parameters  of  control  gate  two  and  three  dc-biases  were 
changed.  This  result  was  expected  based  on  the  DC  observa¬ 
tions  . 

D.  OUTPUT  CIRCUIT  ANALYSIS 

The  operation  of  the  output  circuit  has  been  previously 
described  in  an  earlier  section  of  this  paper.  The  circuit, 
shown  in  Figure  15,  must  provide  an  extremely  high  impedance 
input  to  prevent  loading  the  preceding  stage,  which  in  this 
device  was  the  small  charge  packet  delivered  through  the  chan 
nel  to  the  output  diode.  To  achieve  this  high  input  imped¬ 
ance,  not  only  was  MOSFET  circuitry  employed,  but  the  circuit 
was  capacitively  coupled  to  the  channel  through  the  reverse 


Figure  15 
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biased  output  diode.  While  utilization  of  the  transition 
capacitance  of  a  reverse  biased  diode  does  indeed  represent 
a  high  impedance  coupling  scheme,  it  is  not  without  disad¬ 
vantages.  The  reverse  bias  must  be  maintained.  This  is  an 
additional  requirement  on  the  output  circuit.  The  transi¬ 
tion  capacitance  is  nonlinear;  i.e.,  it  varies  with  the 
applied  voltage.  It  is  the  purpose  of  this  section  of  the 
study  to  examine  these  effects  on  the  output  voltage  as  a 
function  of  the  channel  potential  immediately  adjacent  to 
the  output  diode  within  the  device . 

1 .  Data  Curve  Fitting 

It  became  apparent  that  much  of  the  output  circuit 
analysis  would  require  numerical  solutions  using  a  computer. 
The  output  of  such  techniques  can  be  presented  in  tabular  or 
graphical  form,  but  not  in  closed  analytical  form.  It  was 
desirable  to  have  an  ability  to  reduce  such  data  to  an  ana¬ 
lytical  expression  that  would  closely  approximate  the  orig¬ 
inal  numerical  solution. 

Generally,  a  curve  can  be  closely  approximated  by  a  poly¬ 
nomial.  Such  a  polynomial  would  have  a  set  of  coefficients 
and  be  of  the  form: 

y  =  Cj  +  C2X  +  C3X2  +•  •  •+  C^x*1  1  (8) 

At  the  it^  data  point,  where  the  independent  variable  is  x^ 
and  the  corresponding  dependent  variable  value  is  y^;  the 
error  between  the  polynomial  and  the  measured  data  value  can 
be  mathematically  expressed  as  follows : 
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The  intent  is  to  evaluate  the  coefficients  in  such  a  man¬ 


ner  as  to  minimize  the  total  of  the  errors  at  all  data 
points.  It  is  possible,  in  fact  probable,  that  individual 
errors  could  be  either  positive  or  negative.  In  such  an 
event,  building  a  quantity  based  on  the  sum  of  the  errors 
would  produce  an  invalid  conclusion,  as  large  negative  errors 
added  to  large  positive  errors  could  produce  a  small  total. 

To  generate  a  monotonically  increasing  function  that  is  an 
inverse  measure  of  the  "goodness"  of  fit,  the  sum  of  the 
square  of  the  errors  is  used.  That  is: 


(10) 


(ID 


or, 


To  determine  the  coefficients,  it  is  necessary  to  deter¬ 
mine  the  partial  derivatives  of  the  sum  of  the  square  of  the 
errors  with  respect  to  each  of  the  n  coefficients.  Setting 
these  partial  derivatives  equal  to  zero  yields  a  set  of  n 
simultaneous  linear  equations.  The  solution  vector  of  this 
set  of  equations  is  the  set  of  n  coefficients  what  provides 
a  polynomial  of  degree  n-1  that  best  fits  the  data  points. 


That  is,  the  coefficients  have  been  so  selected  as  to  mini¬ 
mize  the  summation  of  the  square  of  errors  at  each  data  point, 
This  set  of  equations  is  made  up  of  several  summation  terms 
formed  from  the  data  points.  In  the  makeup  of  the  equations, 
wherever  the  summation  symbol  appears,  it  indicates  the  sum¬ 
mation  of  the  variable  term  following  it  for  all  data  points 
from  1  to  m,  the  total  number  of  data  points.  The  set  of 


equations  to  be  solved  is  then: 


C^m 

ClZxi 


+  C0EX.  +  C-EX?  +•••+  C  EX1?"1 
2  i  3  i  n  l 


+  C2ex^  +• 


n 


■+  CnEX1 


Cl£xi 


C.EX1?"1  +  C0EX*?  +« 
1  i  2  i 


-  £y. 

=  ^yixi 
2 

=  Zy.x. 


(12) 


To  facilitate  the  output  circuit  analysis,  a  computer 
program  was  written  to  generate  the  above  set  of  equations 
and  solve  for  the  coefficients  for  any  value  of  n  from  2  to 
10,  with  m  data  points.  For  additional  utility,  the  program 
included  the  capability  to  fit  data  to  a  power  function: 

b 


y  *  ax 


(13) 


and  the  exponential  function: 


bx 


y  =  ae'  (14) 

Both  of  these  functions  can  be  manipulated  for  best  data  fit 
with  the  same  polynomial  approach,  with  the  following  modi¬ 
fications  : 


1.  Use  n  =  2,  that  is  the  best  fit  straight  line, 
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2.  For  the  power  function,  supply  the  natural 
logarithm  of  x,  rather  than  x,  as  the  data, 

3.  Supply  the  natural  logarithm  of  y  as  data, 

4.  Evaluate  the  constant  a,  as  the  exp(C-,)  in 
both  cases,  and 

5.  When  printing  out  the  data,  convert  back  to 
the  original  data  by  using  exp(y.)  and  in 
the  case  of  the  power  function,  exp(x^). 

The  computer  program  that  was  developed  is  appended  here¬ 
to  as  Appendix  D.  In  addition  to  the  form  of  the  curve,  the 
program  will  supply,  as  part  of  the  output,  the  calculated 
values  of  the  dependent  variable,  the  data  points,  and  the 
sum  of  the  square  of  the  errors.  In  a  case  where  no  inform¬ 
ation  is  available  as  to  the  type  of  curve  the  data  repre¬ 
sents,  the  program  is  run  to  try  the  power  function,  the 
exponential  function,  and  all  polynomials  to  degree  9.  The 
sum  of  the  square  of  the  errors  will  be  smallest  for  the  best 
fit  curve . 

2 .  Output  MOSFET  Characteristics 

To  analyze  the  output  of  the  circuit  in  Figure  15, 
a  knowledge  of  the  characteristics  of  the  output  MOSFET  is 
required.  As  previously  discussed  in  section  IV. F. ,  the 
manufacturer  provided  data  on  the  gate  to  source  voltage  and 
the  drain  to  source  current  of  four  sample  devices.  This 
data  is  presented  in  Table  I. 

The  data  for  device  number  3,  was  rejected  as  being  non¬ 
typical  of  the  devices  studied.  The  remaining  data  was  pro¬ 
cessed  using  the  data  fit  routine  described  in  the  preceding 
section.  The  data  was  processed  as: 
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(15) 


yi  *  "W 


xi  =  VDS.  =  VGSi 


n  =  2 


(16) 

(17) 


The  results  of  the  program  are  tabulated  in  Appendix  A 
hereto  and  are  displayed  graphically  in  Figure  10.  The  aver¬ 
age  error  in  the  fit  was  1.29x10”^  amps.  The  equation 
describing  the  current -voltage  relationship  is: 

!DS  =  "k^VGS  "  V 

where:  k  =  2.2032xl0'6  amp/volt2 

Vj  =  -1 .76466  volts 

Recalling  that  the  length  of  the  channel  in  the  direction 
of  current  flow  is  smaller  than  the  designed  or  masked  length, 
L,  because  of  the  p"*"- diffusions  proceeding  laterally  under 
the  mask  by  a  distance  called  the  junction  depth,  x^  .  Then: 


%  =  L  -  2x .  (21) 

J 

In  the  case  of  the  experimental  device: 

W  =  0.4  mils  (22) 

L  =  0.6  mils  (23) 

x.  =  0.06  mils  (24) 

J 

l  =  0.48  mils  (25) 


With  this  geometry,  the  conduction  factor,  k' ,  is: 

k'  =  =  rp  =  2.2644xl0'6  amps/volt2  (2 

^zox 

Using  the  parameters  obtained,  it  was  possible  to  cons- 


(18) 

(19) 

(20) 
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Table  I 


Drain  to  Source  Drain  to  Source  Voltage  -  Vnc. 

Current  -  IDS  VDg  =  VQS  (volts)DS 


(pA) 

Device  1 

Device  2 

Device 

3 

Device  4 

0.1 

2.04 

2.04 

1.84 

1.96 

1.0 

2.50 

2.50 

2.30 

2.44 

4.0 

3.14 

3.14 

2.91 

3.08 

9.0 

3.76 

3.77 

3.52 

3.75 

16.0 

4.42 

4.41 

4.15 

4.37 

25.0 

5.08 

5.10 

4.79 

5.07 

36.0 

5.80 

5.80 

5.47 

5.07 

49.0 

6.46 

6.46 

4.57 

6.50 

64.0 

7.15 

7.16 

6.82 

7.20 

81.0 

7.85 

7.87 

7.49 

7.91 

truct  the  drain 

characteristics  of 

the  output 

MOSFET.  The 

basic  Sah  model 

[Ref.  7]  of 

the  transistor, 

assuming  no 

source  to  drain 

resistance , 

indicates  that  i 

the 

drain  to 

source  current  will  be  a  constant  for  all  values  of  the  drain 
to  source  voltage  once  the  transistor  is  in  saturation.  This 
model  also  provided  expressions  for  the  current  under  condi¬ 
tions  of  saturation,  equation  (4),  and  non-saturation, 
equation  (5) .  The  saturation  relationship  has  already  been 
used  to  define  the  threshold  voltage  and  the  gain  parameter. 
Using  these  values,  the  two  voltage-current  relationships 
were  used  to  establish  the  drain  characteristic  curves  shown 
in  Figure  16.  A  load  line  for  operation  into  a  50,000  ohm 
load  resistance  is  also  plotted  on  the  curves.  This  value 
or  load  resistance  was  typical  of  the  conditions  used  during 
the  experimental  portion  of  this  study. 

3 .  Output  Diode  Effects 

Capacitively  coupling  the  output  MOSFET  to  the  CCD 
channel  through  a  reverse  biased  diode  provides  the  required 
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Figure  16 
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high  impedance  load  for  the  channel.  Unfortunately,  it  pre¬ 
sents  the  complication  of  having  the  capacitance  of  the  coup¬ 
ling  capacitor,  the  diode,  vary  as  a  function  of  the  voltage 
across  it.  Asstiming  that  the  reverse  biased  junctions  invol¬ 
ved  with  the  output  circuit  are  made  with  abrupt  changes  from 
n-region  to  p-region,  the  value  of  the  junction  capacitance 
is  [Ref.  8] : 


2 


(27) 


farad/cm 


-19 

where:  qQ  =  electronic  charge,  1.602x10  coulomb 


er  =  relative  permittivity  of  the  oxide,  11.8 

-1 4 

eQ  =  permittivity  of  free  space,  8.854x10  farad/cm 

3 

=  concentration  of  donors  in  n-region,  atoms/cm 


V  =  applied  potential  in  volts. 

a 

To  assist  in  analyzing  the  effect  of  the  output  diode, 
the  simplified  circuit  representation  in  Figure  17  was  used. 
The  input  potential,  that  potential  adjacent  to  the  cathode 
of  the  output  diode,  is  the  signal,  x.  The  current  i-^,  flows 
through  the  output  diode  and  then  branches  into  currents  i ^ > 
i^  and  i^.  Current  flows  through  the  reverse  biased  p-n 

junction  represented  by  the  reset  MOSFET,  Q5 .  Current  i^, 
represents  the  leakage  through  the  substrate,  the  input  resis¬ 
tance  of  the  output  MOSFET  and  the  pinched  off  channel  resis¬ 
tance  of  the  reset  MOSFET.  Current  i^,  represents  the  current 
into  the  gate  of  the  output  MOSFET,  Q6.  The  signal  output  is 
,  the  input  voltage  to  the  output  MOSFET.  Writing  the  nodal 
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V^T 


EQUIVALENT  CIRCUIT  SEEN  BY  THE  OUTPUT  DIODE 

Figure  17 


equations  for  node  A: 


Now,  in  general 


il  =  '2  +  ’3  +  ’4 


(28) 


'  =  5t<cv> 

Then,  for  the  circuit  in  Figure  17: 

k,  ,  k,  dV.  k,  dV.  V.  dV, 

ldx  1  l  2  1  l  ^  r 

2;X-Vi  Ht  "  2/x-V.  dt  "  2/v.-Edd  ^t"  Rt  LTdt 


(29) 


(30) 


or,  solving  for  the  derivative  of  the  output  with  respect  to 
time :  ,, 


1  dx  1 


dt 


2 A  -  T.  dt  Y'i 


+  CT  + 


(31) 


where : 


=  Ai/ 


9o£r£oNd 


(32) 
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(33) 


.  A 


CT  = 


e  e  A- 
r  o  3 


(34) 


Some  of  these  quantities  are  known,  but  for  others  assump¬ 
tions  must  be  made.  Every  assumption  presents  a  possible 
error  in  the  analysis  of  the  experimental  circuit,  but  the 
values  assumed  will  be  typical  of  those  seen  in  semiconductor 
technology.  Further,  if  the  assumptions  depart  from  the 
experimental  device,  they  will  still  provide  meaningful  ana¬ 
lysis  of  the  type  of  output  circuit  employed. 

The  value  of  the  oxide  thickness,  d  in  equation  (34)  is 
known  as  1000A.  Likewise,  the  values  of  qQ,  N^,  er  and 

e ^  are  established  quantities.  The  initial  conditions,  the 
input  voltage  function,  x,  and  the  areas  A^,  A ^  and  A^,  are 
not  so  well  defined. 

The  area  of  the  output  diode,  A-^,  is  taken  to  be  a  rec¬ 
tangle  2.0  mils  long  by  0.25  mils  wide.  The  area  of  the  reset 
MOSFET  source  diffusion  is  taken  to  be  0.4  mils  long  by  0.25 

mils  wide.  The  output  MOSFET  gate  is  0.4  mils  by  0.6  mils. 

9 

The  distributed  resistance,  R^.,  is  assumed  to  be  5x10  ohms. 

Since  the  quiescent  output  voltage  was  experimentally  seen 
to  be  about  -12  volts,  it  can  be  seen  from  Figure  16,  that  the 
quiescent  gate  voltage  must  be  about  -24  volts.  The  initial 
value  of  the  variable,  V.  ,  is  then  taken  to  be  -24.0  volts. 

i 

The  input  function  is  taken  to  be  a  voltage  step  from 
-25  volts  to  -12  volts  at  time  equal  15  usee  returning  to 
-25  volts  at  the  time,  47  usee. 
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Using  the  known  values  and  the  assumed  quantities,  the 


numerical 

constants  now  become: 

k.j  =  2.2951x10"^  farad-volt^ 

(35) 

k2  =  5. 902x1 0-1 5  farad-volt^ 

(36) 

C-j.  =  1.618x10”^  farad 

(37) 

A  computer 

program  to  solve  the  general  initial 

value  prob- 

lem: 

at  ■  f(x'  y-  z> 

(38) 

where : 

x  =  g(t) 

(39) 

z  = 

dt 

(40) 

and : 

<< 

11 

<< 

O 

(41) 

when : 

t  =  0 

(42) 

was  developed  using  fourth  order  Runge-Kutta  methods.  This 
program  is  listed  in  Appendix  B  hereto.  The  program  was  exe¬ 
cuted  and  the  results  are  tabulated  in  Appendix  A  and  graph¬ 
ically  displayed  in  Figure  18. 

The  results  of  the  solution  to  the  differential  equation 
represent  the  response  of  the  output  diode  to  an  input  pulse 
in  the  CCD  channel.  This  response  would  represent  one  sample 
of  the  output  signal.  In  order  to  determine  the  linearity  of 
the  output  diode  to  an  input  waveform,  of  which  the  function 
g(t)  is  only  one  sample;  and  which  will  be  processed  by  a 
sample  and  hold  circuit,  it  was  necessary  to  modify  the 
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DIODE  OUTPUT  VOLTAGE,  v,  (VOLTS) 

-24.0  -218  -23.6  -23.4  -23,2  -23.0  -22.8  -22.6  -22.4 


Figure  18 
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-30  -25  -20  -15  -10 

DIODE  INPUT  VOLTAGE,  x  (VOLTS) 


computer  program. 

This  was  done,  and  the  results  are  discussed  as  the  input 
to  the  source  follower. 

4 .  Input  to  the  Source  Follower 

While  the  transient  analysis  of  the  output  diode  to 
a  step  function  provides  the  waveform  of  the  input  to  the  out 
put  MOSFET,  it  does  not  represent  the  response  to  the  signal 
input.  Rather,  it  represents  the  response  to  a  single  samp¬ 
ling  of  the  signal.  What  is  required  is  the  relationship 
between  the  peak  amplitude  of  the  potential  adjacent  to  the 
output  diode  in  the  CCD  channel  and  the  peak  of  the  output 
diode  response. 

To  accomplish  this,  the  differential  equation  solution 
routine  was  modified.  The  program  was  modified  to  solve 
equation  (31)  forty-one  times.  The  amplitude  of  the  chan¬ 
nel  potential  step  function  was  increased  from  -20  volts  to 
-10  volts  in  0.25  volt  steps.  The  duration  of  the  step  and 
the  quiescent  value  of  the  voltages  were  maintained  as  before 
The  program  was  required  to  provide  an  output  of  only  the 
peak  input  pulse  amplitude,  x,  and  the  peak  output  pulse, 

V.  . 

i 

The  modified  program  is  appended  hereto  in  Appendix  C. 

The  results  of  the  program  are  tabulated  in  Appendix  A  and 
presented  graphically  in  Figure  19.  The  solution  data  from 
this  program  was  used  as  an  input  to  the  curve  fitting  rout¬ 
ine.  It  was  found  that  the  best  fit  of  the  results  was  with 
a  polynomial  of  degree  3.  The  expression  for  the  peak  MOSFET 
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PEAK  DIODE  OUTPUT  VOLTAGE,  Vj  -  (VOLTS) 
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Figure  19 
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input  voltage  as  a  function  of  the  peak  channel  potential 
adjacent  to  the  output  diode  is: 

-V.  =  22.15  +  3.81xl0-2tx|  -  6.15xl0'4|x|2  +  2.33xl0"5|x|3  (46) 

The  average  error  of  fit  was  found  to  be  2.41x10  ^  volts. 

An  advantage  of  presenting  the  results  as  a  polynomial  is 
that  a  measurement  of  linearity  is  available.  The  value  of 
the  coefficient  of  x  to  the  first  power  may  be  ratioed  to  the 
sum  of  the  coefficients  of  all  powers  of  x,  to  express  the 
fraction  of  linearity.  Likewise,  the  sum  of  the  coefficients 
of  all  x  terms  with  exponents  larger  than  one  is  an  expres¬ 
sion  of  the  departure  from  linearity.  Using  the  more  precise 
coefficient  values  tabulated  in  Appendix  A,  it  is  seen  that 
the  output  of  the  output  diode  is  98.35  percent  linear.  Or, 
only  1.65  percent  of  the  output  variation  is  caused  by  non¬ 
linear  relationships  to  the  input. 

5 .  Output  of  the  Source  Follower 

The  next  element  in  the  output  circuit  is  the  output 
MOSFET  itself.  Using  the  model  of  the  transistor  that  has 
been  used  throughout,  an  expression  for  the  output  as  a  func¬ 
tion  of  the  MOSFET  input  can  be  developed.  Referring  to 
Figure  20  for  nomenclature,  the  expression  for  the  output  cur¬ 
rent  as  a  function  of  the  gate  input  voltage  will  be  developed. 
By  examination  of  Figure  16  remembering  that  the  output  quies¬ 
cent  voltage  was  about  -12  volts,  it  is  seen  that  with  a  drain 
bias  supply  of  -25  volts,  that  the  transistor  is  always  in 
saturation.  Utilizing  equation  (18),  the  output  current  is 

seen  to  be : : 
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(47) 


but , 


then , 


(49) 


(48) 


and  by  simple  manipulation, 


V0  -  -  !<2<Vt  -  vf)  ♦  i^)±  Kt2<VT-Vi)+W,i!-4(V^'4-2V,»T) 


(50) 


This  output  relationship  was  evaluated  using  the  computer 
with  the  program  listed  in  Appendix  E.  The  input  voltage,  V^, 
was  programed  into  the  calculation  as  the  third  degree  poly¬ 
nomial  developed  in  the  preceding  section.  The  negative  opt¬ 
ion  for  the  radical  in  equation  (50)  was  ignored  in  the  eval¬ 
uation,  since  it  would  produce  output  voltages  more  negative 
than  the  drain  bias  voltage,  E^,  a  physical  impossibility. 

The  results  of  the  computer  run,  then,  was  a  tabulated 
listing  of  the  output  voltage  at  the  source  follower  load 
resistor  as  a  function  of  the  channel  potential  adjacent  to 
the  output  diode  within  the  CCD.  This  is  the  relationship 
that  was  desired  as  a  result  of  the  output  circuit  analysis. 
The  results  are  tabulated  in  Appendix  A,  and  presented  graph¬ 
ically  in  Figure  21.  Visually,  the  output  voltage  dependence 
on  the  channel  peak  potential  appears  quite  close  to  a  linear 
relationship . 

Once  again,  the  curve  fitting  program  was  employed  to 
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SIMPLIFED  OUTPUT  CIRCUIT 
Figure  20 

determine  an  analytical  expression  for  the  output  voltage 
pulse  peak  as  a  function  of  the  peak  potential  in  the  CCD 
channel  adjacent  to  the  output  diode.  The  best  fit  was  found 
with  a  polynomial  of  degree  3.  The  functional  relationship 
is  : 

-V  =  10.58  +  2 . 72x1 0-2lxl  -  4.98xlO"V  +  1 .78x10"5ixi3  (51) 

0 

The  average  error  was  found  to  be  6.77x10  volts. 

Once  again  using  the  more  precise  coefficients  in  Appendix 
A,  the  output  circuit  is  98.14  percent  linear.  That  is,  only 
1.86  percent  of  the  variation  in  the  peak  output  voltage  due 
to  the  channel  potential  pulse  peak  is  contributed  by  higher 
order,  or  nonlinear,  terms. 
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PEAK  SYSTEM  OUTPUT  VOLTAGE,  V0  -  (VOLTS) 


PEAK  CHANNEL  POTENTIAL 
ADJACENT  TO  THE  OUTPUT  DIODE ,  X  -  (VOLTS) 

Figure  2  1 
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E .  LINEAR  MODEL 


If  electronic  devices  can  be  operated  so  that  the  non¬ 
linear  distortion  is  small,  then  they  can  often  be  approxi¬ 
mated  by  linear  equations  and  circuits.  In  such  cases,  all 
the  procedures  of  linear-circuit  analysis  can  be  applied  and 
complex  circuits  can  be  analyzed  with  a  minimum  of  effort. 

In  this  section,  a  linear  model  of  the  CCD  in  an  analog  sig¬ 
nal  processing  mode  is  suggested,  but  no  experimentally 
determined  parameter  values  are  provided. 

The  currents  of  the  network  in  Figure  22  can  be  repre¬ 
sented  as  linear  functions  of  the  voltages.  Thus, 


=  ynvi  +  yi2v2 

(52) 

=  y21vi  +  y22V2 

(53) 

The  parameters  of  these  equations  are  termed  y  parameters. 
They  have  the  dimensions  of  admittances.  It  is  most  conven¬ 
ient  to  measure  these  parameters  under  short-circuit  condi¬ 
tions.  It  was  the  experimental  application  of  this  concept 
that  proved  beyond  the  capabilities  of  this  experimenter. 

The  intent,  was  to  apply  a  voltage  generator  to  the  input 
gate  and  short-circuit  the  output,  or  vice  versa.  Or,  fail¬ 
ing  that,  with  the  circuit  otherwise  operating,  to  determine 
the  parameters  from  the  definitions,  that  is: 

aIl 

yu  =  jy1’  Short-circuit  input  driving  point  admittance  (54) 


'1 
3  I 

_  _ 2  Short-circuit  forward  transfer  admittance 

y21  “  3V-,’ 


(55) 
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(56) 


y12  = 

3 1 2 

y22  = 


» 


5 


Short-circuit  reverse  transfer  admittance 
Short-circuit  output  driving  point  admittance 


(57) 


In  working  with  the  CCD,  it  became  apparent  that  little 
or  no  short-circuit  reverse  transfer  admittance  existed.  In 
light  of  the  manner  in  which  the  CCD  is  fabricated  and  func¬ 
tions,  this  is  not  unreasonable.  It  was  also  apparent  that 
the  short-circuit  input  driving  point  admittance  was  small 
and  mainly  capacitive.  The  arrangement  of  the  output  circuit 
would  indicate  that  the  short-circuit  output  driving  point 
admittance  would  be  the  reciprocal  of  the  output  MOSFET's 
dynamic  drain-source  resistance.  The  quantity  that  would  be 
peculiar  to  a  CCD  was  the  short-circuit  forward  transfer 
admittance,  Y2i-  This  quantity  would  have  to  indicate  the 
sampled  nature  of  the  output,  as  well  as  the  time  delay  of 
the  output  with  respect  to  the  input.  Neglecting  the  output 
dynamic  drain-source  resistance,  such  a  quantity  would  be: 


y21 


dl 


dVi 


u(t  -  5^-)  -  u(t  -  ir-  -  x) 


(58) 


where : 


dl 

=  The  slope  of  the  dc  gain  curve. 

u(t)  =  The  unit  step  function 

n  =  The  sequence  of  consecutive  integers  that 
best  represents  the  sampling  during  the 
period  of  interest. 


N  =  The  number  of  charge-coupled  structures 
f  =  The  phase  clocking  frequency 
t  =  The  duration  of  the  sampling  pulse. 
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In  arriving  at  N,  the  number  of  charge-coupled  structures,  a 
structure  constitutes  one  delay  element,  that  is,  in  a  two- 
phase  device,  a  set  of  four  consecutive  gates  connected  to 
the  two  different  phase  clocking  signals. 


In  those  applications  where  multiple  inputs  are  used,  a 
modification  would  have  to  be  imposed  on  the  model.  Each 
input  would  have  its  own  input  driving  point  admittance;  but 
beyond  that  point,  each  input  would  have  to  be  ratioed  based 
on  its  effect  on  the  output  current.  The  ratioed  potentials 
would  act  on  a  summing  point,  the  voltage  of  which  would  be 
the  input  voltage  to  the  expression  for  the  model's  dependent 
current  source. 
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VI.  CONCLUSIONS 


It  was  established  that  the  two-phase  charge-coupled 
shift  register  could  be  used  to  process  analog  signals.  In 
the  course  of  such  processing,  a  consistent  gain  was  dis¬ 
played,  and  an  electronically  controlled  delay  could  be  insert¬ 
ed.  The  output  of  the  chrage-coupled  shift  register  was  a 
sampled  version  of  the  input  signal,  but  a  sample  and  hold 
circuit  restored  the  wholly  analog  nature  to  the  output  sig¬ 
nal.  Little  nonlinear  distortion  was  observed  over  a  small 
range  of  input  signal  level,  i.e.,  0.1  to  0.2  volts. 

Further,  it  was  found  that  the  charge-coupled  shift 
register  could  be  made  to  function  as  an  analog  signal  adder- 
sub  tractor  with  good  linearity.  Once  again  the  dynamic  range 
over  which  linearity  could  be  achieved  was  0.1  to  0.2  volts. 
The  small  dynamic  range  is  attributed  to  the  small  surface 
potential  wells  of  the  experimental  device.  If  more  charge 
could  be  transferred,  a  larger  dynamic  range  would  be  evi¬ 
denced. 

An  analysis  of  the  output  circuit  employed  established 
that  very  slight  nonlinearities  were  introduced  by  the  out¬ 
put  technique  used  on  the  experimental  device. 

Finally,  a  linear  model  of  the  charge-coupled  shift 
register  is  proposed. 
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TABULATION  OF  COMPUTER  OUTPUTS 
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-  5  *d 

-  4  .  <4  4  b 

-5*9  69 72  E -3 

-  o  •  4  h 

-  #447 

-6.96937E-3 

-  0  •  4  h 

-  4  7 

- b  •  9  b 4 3  7 E-3 

-  h  •  5 

-  4  .42  7 

-  7  •  4  2  o  7  4  F  -  3 

-7.1b 

-4  .V*  4  8 

-  7.993  54F-3 

-  7  •  1  o 

-4  •  4  '2  d 

-  d  .  4  /)  6  3  9  E  -  3 

-7.2 

-4  . 2(4 b 

-d *4F  775E-3 

-  7  *d  5 

--2  •  /  4  9 

-9.432  5  7E-3 

-  7  .d  7 

-4  *449 

-9  *4b22  5E-3 

-  7.3  1 

—  2  •  4  '4  9 

-9.12  1  F3E-3 
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i.  c 

•  J  a4 

J  •^dd4dt-b 
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j  •  4  9  -?  4  9  t  -  o 

-12 

.  J  "v 
f.  c. 

•  9  a  2  3 
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D  I  DL01  A 


2  2  :  1  b 


03-DEC-75 


p£AK  OUTPUT  DluOE  VOLTAGE  AS  A  FUNCTION  OF 

peak  channel  voltage 


X  PEAK 

Y  PEAK 

-10 

-22  .4928 

-10.25 

-22 .501 

-Id. 5 

-22.5091 

-Id. 75 

-22*51  73 

-11 

-22 .5256 

-1  1  .25 

-22.5338 

-11*5 

-22  .542  1 

-11.75 

-22  »5  5d  4 

-12 

-22  .5563 

-12.25 

-22 .5o72 

-12.5 

-22.575c 

-12-75 

-22  .  58  4 

-13 

-22 .5925 

-13.25 

-22 .60  1 

-13.5 

-22  «  6  d  9  6 

-13 .75 

-22  .61o2 

-14 

-22  .6263 

-1  4.25 

-22  *63  55 

-14.5 

-22.6443 

-14.75 

-22  .653 

-15 

-22  *661 o 

-1  5*25 

-22 .6  70  7 

-15.5 

-22.6  796 

-15.75 

-22  .638  b 

-16 

-22 .69  76 

-  1  b  *2  5 

-22 • 7d  6  7 

- 1  b  *  5 

-22 .71  5a 

-16.75 

-22  •  725 

-1  7 

-22.73  43 

-1  7.25 

-22  .  743  6 

-17.5 

-22 .753 

-17.75 

-22  •  762  5 

-18 

-22 .7721 

-18  .25 

-22 .78  18 

-la  .5 

-22  *79  1  6 

-18.75 

-22  .30  1  5 

-Id 

-22  .3115 

-18*25 

-22  .32  1  6 

-19.5 

-2  2  .33  19 

-19.75 

-22  .8  423 

-2d 

-22  .3  529 

TIME:  Idol. 21  SECS. 
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PEAK  DIODE  OUTPUT  VOLTAGE 
As  a  Function  of 

PEAK  CHANNEL  SIGNAL  POTENTIAL  ADJACENT  TO  THE  DIODE 


K  U  i\ 

(Least 

Squares  Curve  Fit) 

ol  bU>2 

£  V1  2  4  4 

0  4 -DEC- 7b 

Qur^UT  VOLTAGE 

.*5 

a  flmc  r  i 

00  u  r 

1  i\r  {J  T  VOLTAGE 

H  OL  r  .tO.il  AL  CUk'VE 

r  1  f : 

Y=(  22*15  )♦( 

3.9394, )E-2  ) 

»  A  ♦ 

(-b»1olr?5c-4  ) 

*( 

2  * 

33  1  33E-o 

)  *a  t3 

♦( 

/) 

)  *  A  T  5  (  t) 

)  *  X  T  •  •  • 

♦( 

/) 

) *  X  f  7+ (  H 

)  *  X  f  8  ♦  •  •  • 

♦( 

6 

)  *  x  r  3 

DATA  PCI  NTS  A riD  CALCuLAfcO 

VALUES  OF  Y 

VALUE  OF  A 

VALUE  U^  Y 

C  ALC.JLA  I  EG  Y 

22 

•  -13  2  o 

22 

.432  7 

1 6  »  2  b 

22 

.  5/1 

22 

•  5/ Hi 

12*5 

22 

•  5/131 

22 

•  5u  i2 

1/1.75 

22 

.51  73 

22 

.51  74 

1  1 

22 

*52  56 

22 

•  52  5o 

1  1  .2  5 

22 

•  53  38 

22 

*5339 

1  1  .5 

22 

.  542  1 

22 

.  5422 

11.75 

22 

*55/14 

22 

» 5  5 .1  5 

12 

22 

.  55o8 

22 

•  55d9 

12.25 

22 

*  5b  72 

22 

•  5  o  72 

12  • 5 

ic. 

*  5  7  bo 

22 

•  5  75b 

12-75 

22 

*59  4 

22 

•  58  4 

13 

22 

•  592  5 

2  2 

.  5  3  2  5 

13*25 

22 

•  oil 

22 

•  «/}  1 
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Value  of  X 

Value  of  Y 

Calculated 

1  a  .  5 

2  2  *  n  t)  9  o 

22  •  6<>S  5 

1  3  *  7  o 

2  2  » o  1  o  2 

2  2  *  6 1  a  1 

1  4 

22  * 52  t>5 

2  2  *62  8  7 

14*25 

22 .6355 

22.63  54 

14.5 

2  2 . 6  44 J 

22.6441 

14.75 

2  2  *  6  53 

22 .6529 

1  5 

2  2  *  h  o  1  b 

2  2  » 6  b  1  7 

15.2b 

2  2  •  -  7 ,  7 

2  2  » 6  7  v;  0 

1  3 » 5 

22  *o  7jb 

22  *  8  7H  5 

1  o  *  7  5 

2  2  •  H  o 

2  2  *  ^  3  3  5 

1  0 

2  2  .E-?  7  b 

2  2  »  6  ?  7  3 

1  o .  2  5 

22  *  7  5  o  7 

2  2  .  7  2'  r  c 

1  c  *  5 

22 . 71  bo 

2  2  *  7 1  3  a 

1  o  *  7  5 

2a  .  725 

22  *  72  o  1 

1  7 

2a  •  73  4J 

22 .73  44 

17.25 

22  •  74 j  0 

22.7 43a 

17*5 

22 • 752 

22  •  7  5  J  2 

17.75 

2  2  »  7b2  5 

22  •  7c2 7 

1  0 

22  .  7  72  1 

22.7/2  4 

1  3  »2  3 

22  »  73  1  0 

22 .73  2  1 

19*5 

22  •  7  3 1 b 

2  2  *  7  i  1 0 

1  <3  *  7  5 

2  2  *  9  a  1  5 

22*3217 

19 

22.6115 

22  .9  1  1  0 

15*25 

22  .82  1  c 

22  .32  1  7 

19.5 

22  .83  19 

22  .9  3  13 

19.75 

22  .9  423 

22  .9  42 

2  4 

22  .9529 

22  .9  523 

THE  a  IJ'*  Or  TrtP  3  0U48E  OF  f  r  E  Err.  uFS  -  n. 
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OUTPUT  OP  SOURCE  FOLLOWER  OS  0  FUNCTION  OF 
CHANNEL  VOLTAGE  IN  CCO 


C  rl  A  N  N  E  L 

OUTPUT 

1  N 

VOLTAGE 

VOLTAGE 

SATURATION 

-1/4 

-10*3182 

YES 

-10  -2 5 

-10  .3239 

Y  ES 

-10.5 

-10.829  5 

YES 

- 1 o  *  75 

-10  .33  52 

YES 

- 1 1 

-10  .3  403 

YES 

-1  1  .25 

-10*9465 

YES 

-11.5 

-10  .3522 

YES 

-11.75 

-  10  .3  5  79 

Y  ES 

-12 

-10.363b 

Y  Eo 

-12  .25 

-  1  0  *3  b3  3 

Y  ES 

-12.5 

- 1  0  *8  7  51 

Y  cS 

-12*75 

- 1  0  .fldoy 

YES 

-1  J 

-  1  0  *8  8  b  7 

YES 

-13.2  5 

-  10  .892  5 

Y  ES 

-13.5 

-  1 0  *3  3o  4 

Y  Ed 

-13.75 

-10 .OU  43 

YES 

-1  4 

-10  .9  1/2 

YES 

-14.25 

-10  .91  o2 

YES 

- 1  4 . 5 

-10  .92  2  1 

YES 

-14.75 

-10 » 9  2  o  2 

YES 

-15 

-10 .9342 

Y  ES 

-15.25 

-10  *3403 

Y  Ed 

-15.5 

-  10  *  9  4  b  5 

YES 

-15.75 

-10  .952b 

YES 

-Id 

-10  .9  58  8 

Y  ES 

- 1 o»25 

-10  *9651 

YES 

- 1  b  *  5 

-10.9  714 

Y  ES 

in 

r*- 

• 

o 

1 

-  1  0  .9  7  78 

Y  ES 

-17 

-10  .98  42 

YES 

-1  7.25 

-l/i  »99ob 

YES 

-17.5 

-10.-9  71 

YES 

-1 7. 75 

-1  1  .003  7 

Y  ES 

-18 

-1  1  .0  103 

YES 

-13  .25 

-1 1  .0  1  7 

Y  ES 

-  1  d  »  5 

-1  1  .023  7 

Y  ES 

-18.75 

-  1  1  *  0  3  0  5 

Y  Ed 

-13 

-  1  1  .03  73 

Y  ES 

-13.25 

-  1  1  .0  4  42 

YES 

-19.5 

-1  1  .0512 

YEd 

-1y . 75 

-11 *0  5b  2 

YES 

-2  0 

-11  .Oo53 

YES 

T I  1*02  S  ECS • 
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PEAK  SYSTEM  OUTPUT  VOLTAGE 
As  a  Function  of 

PEAK  CHANNEL  POTENTIAL  ADJACENT  TO  THE  OUTPUT  DIODE 
K'Jts  (Least  Squares  Curve  Fit) 


8  1  JL  Cl 

2  2:11 

f>5- DEC -75 

OUTPUT  VOLTAGE 

P  9  u  M  vOSFET 

SGUhCc  FuLL0*E8 

AS  A  plJl'iO  T 1  ON  OF  ChA.\,\E 

L  V  0  L  T  A  G  E  . 

HOLY  i\ 0  'i  1  A  L  C  j- 

V  6  FIT: 

Y  = (  1 /  •  5 7 o  )*( 

2  -  72  1  3  a  £ -2  ) 

•  a  +■  (-4.376588 

1.775966-5  ) 

*  A  T  3  —  (  2  )  *  A  f  4  + 

♦( 

0  )«Xt5>(  2 

)  *  x  f  b-*’  ... 

♦< 

2  )  *  A  f  7*  (  2 

)  *  A  i'  8  +  .  .  . 

*( 

2  )  *  A  1 9 

DA  TA  r 0 1 N T 3  A\C  CALOULA  TEu 

VALUEi  OF  T 

VALUE  OF  X 

value  of  r 

CALCuL A  TED  Y 

1  2 

1 2  .8  182 

12  .3  182 

1  J  -  2  5 

10.3239 

1 2  »3  238 

12.5 

10.8295 

12  .8295 

1  vi  .75 

12.33  52 

12  .33  51 

1  1 

12 .8  42  8 

1 2  .3  42  8 

11.25 

1  2  •  8  4  b  5 

1  2  *3  4 b 5 

11.5 

10  .9  52  2 

1  2  .3  522 

11.75 

10.3  5  79 

1  2  *3  5  78 

12 

12 . 8  6  3  o 

12  *8  b3  b 

12.25 

1 0 .3  b93 

1  ’<?  *3  69  4 

12.5 

12 .3  751 

12  .8  751 

12.75 

12*8509 

12.3829 

12 

12*38  67 

1  2  » 8  8  b  7 

13  .25 

12  .892  5 

12  *3926 
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Value  of  X 

Value  of  Y 

Calculated 

13*5 

10  .396  4 

1  b  « 8  9  6  4 

13.75 

1  v!  .9  0  43 

10*9043 

1  4 

10  .9  102 

10.9102 

1  4.3  0 

10 .9 1 62 

1  0.91  62 

14.5 

10 .922  1 

10*9224 

14.75 

1*'  *  9  2  8  2 

10  .9  2 a 2 

1  5 

10  .93  42 

1  /  .9  3  42 

15.3b 

1  0  .9403 

10.3403 

15.5 

10.2465 

1  0 . 9  4  r  4 

1  5.  75 

1 0  *  9  52  o 

1  0 . 9  o2  6 

1  6 

10.9  563 

1  '0  .4  5o  d 

16.25 

1  0  *9651 

1  0  .  9  D  5 1 

16*5 

10*9  714 

10.9714 

16.75 

1  0  .9  7  73 

1  4  .9  77  7 

1  7 

1  '4  . 9  o  4  2 

10  •  9  n  4  1 

1  7.25 

10.9906 

1  '0  *  4  3  f  i  0 

1  7.5 

1  0.99  71 

10.49  71 

17.75 

1  1  .0037 

1  1  . 0  0  3  o 

16 

1 1  .0  103 

11.4103 

18.25 

1  1  .01  7 

1 1 .01 69 

13.5 

1  1  .023  7 

1 1  .023  7 

13*75 

11.0305 

1  1  .0305 

19 

1  1  .03  73 

1 1 .03 73 

19.25 

1 1 .0442 

11.0442 

19*5 

1  1  .0512 

11.0512 

19.  75 

1  1  .0562 

1  1  .0  563 

2H 

1  1  .0603 

1  1  .0654 

f  H  £  SUM  uF  The  SQUAne  Or-  I  rt  £  E«*Gr<S  =  7  .  70  5  1 
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3  I  CL 1 


21:34 


03 -DEC  -  75 


10  0 ' 

N  A  ^  E  — 

d 1 DL4  1 

110  * 

120  ' 

DE3CP  1 

P  T  1  ON  —  SOLVES  THE  INITIAL  VALUE  PR  03 LE  4  * 

130  ' 

1  40  ' 

Y  '  =  F 

(  <  *  Y »Z) 

1  54  ' 

whEKE"  X  =  G(  T) 

160  ' 

Y ( TU)=Y0 

1  70  ' 

Z( T)=X  ' 

180  ' 

190  ' 

SOURCE 

-  -R  UNoE  -K  U  T  T  A  (FOURTH  ORDER  ACCURACY) 

200  ' 

210' 

1  NS  Thu 

CTI0N5--THE  INTEGRATION  STEP  IS  H  •  AM  THE 

2  2  0' 

SOLUTION  IS  TABULATED  OVER  THE  INTERVAL  T4>=T>=b* 

2  3  0  ' 

IN  S  T  E 

r'S  OF  SIZE  L  • 

2  40  ' 

2  50  ' 

I  £  FOLLOWING  PROGRAM  1  N  r  U  T  5  A  £  R  E 1 U  1  R  E  0 : 

2  ot;  ' 

0  70  ' 

514 

PR  1  NT  "DES IRED  OUTPUT  Tl  TLE" 

200  ' 

530 

DEF  FNF( x * Y *Z) =  0**** ] 

29  0  ' 

610 

DA  TA  T4  *  Y4  *  6 

3  0  0  ' 

620 

0  A  T  A  rl  »  L 

314  ' 

640 

LET  m=0.0  OH  1.0 

324  ' 

09  0 

0  E  r  FNG( T ) =C  *****  3 

33w  ' 

3  40'  THE  OUTPUT  *111  3  E  THE  TABULATION  OF  T*  X*  AND  Y  t-RGtf 
3  54*  f  0  10  3*  IF  4  =  0  •  2  •  IF  ,8=1*0  IS  SELECTED*  THE  VALUES 

35k)'  OF  X'  AND  V  -ILL  ALSO  3E  TABULATED. 

3  70  ' 

380  '  FOR  MODIFICATIONS  OF  HE  PRINT  ROUTINE*  F  OH  EXAMPLE  TO  PLOT 
390  '  I H  E  SOLUTION  HA  HE*  THAN  TO  TABULATE  ITS  VALUES*  HE  FCLLO  a  I  N 

40  a  '  LINES*  AT  LEAST*  SHOULD  3  E  “GOIFIED: 

4U  ' 

420'  850  Pul  NT  “VALUE  OF  T''*"  VALUE  OF  X“*“VALUE  Or  Y" 

43  0'  8  60  PRINT 

440'  900  PRINT  T*  X*  Y 

450'  1  4  bU  HR  I  NT  T*  X3*  Y7 

4  6.)  ' 

470  ' 

450 '  *******  MAIN  PROGRAM  ********* 

490  ' 

500  REM  INSERT  THE  PRINTOUT  TITLE* 

510  PRINT  "CCD  OUTPUT  DIODE  RESPONSE  TO  A  STEP  FUNCTION  INPUT" 

511  PRINT 

520  REM  DEFINE  HE  DIFFERENTIAL  EOJATION. 

530  DtF  F N F ( X  * Y  * Z  ) 

540  LET  C=1.4755£-14/SQR(APS(x-Y)) 

550  LET  D=2.951E-15/SQR(A8S(Y>25.0)  ) 

5o0  LET  F=C*Z-(2.0E-10)*Y 

b  7'0  LET  G=C* D+  1  •  6 1 8 E~  1  3 

580  FNF  =  _F/G 
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59  2  FNENO 

bv?^  REM  INSERT  ITERATION  AND  TABULATION  DATA. 

6l3  'Data  1  •  a 9 a 5 E - 5 ♦  -24. 9*  4.7915E-5 
62 2  DATA  5.9E-9*  1.0E-6 

639  REM  SPECIFY  THE  TYPE  OE  OUTPUT. 

o  4  9  LET  R  =  0  •  0 

65w  READ  TO*  Y*  6»  H*  L 

669  LET  E  =  1.0/L 

679  LET  T=T0 

690  .  REM  OEFINE  THE  OR  I  VI  NG  FUNCTION 

e90  OEF  FNG( T ) 

799  FNG=-25»3 

719  IF  T<15»9E-6  Th£N  749 

729  IF  T  >4  7  *0  E-6  THEN  749 

73  2  F N G  =  - 1 2  •  0 

7  49  F  N  E  N  0 

759  I f  (d-T2)*H>0.0  THEN  773 

7  69  LET  H  =  -H 

7  79  I  r  (  d  -  T  9  )  *  L  >  9  •  9  T  h  E  N  792 

769  LET  L=-L 

79  3  LET  X  =  F  N  6  (  T ) 

d  3  9  LcT  Zs(FNG(T+-h/2)-FnG(T-h/2))/H 

d  1 9  K  t  M  PRINT  INITIAL  CONDITIONS, 

d IF  M  =  9.0  THEN  9  59 

d  3  3  P  R  I  N  f  "  T  "  * "  <  »,  i.  y  *•  » 11  X  '  ••  ,  >•  y  1 

9  42  GO  T  0  9  62 

d 59  PkI NT  "VALUE  OF  T " » " VALUE  OF  X">"VALUE  OF  Y" 

3  69  PRINT 

879  IF  M  =  9  *  9  ThEN  993 

6d  9  PRINT  Tt  X»  Y  »  2*  rNF(X»Y»Z) 

d99  GO  TO  919 

999  PhINT  T  »  A  »  Y 

919  GO  TO  1999 

922  REM  PERFORM  THE  aASIC  ITERATION. 

939  LET  Z=(FNG( T*H/2)-FNG( T-H/2))/H 
949  LET  x  =  F N G (  T) 

359  LET  K 1 =FNF( X  * Y  »Z ) 

9  b9  LET  Z  1  =(  FN6(  T  +  H  ) -FNG(  T )  ) /H 

9  7 id  LET  X  1  =  F’NG(  T-.H/2  ) 

989  LET  Y1=Y*H*Kl/2 

999  LET  K2=FMF( X  1  *Y  1  *Z1 ) 

1999  LET  Y2=Y-*-H*K2/2 

1919  LET  K3=FNF(X1 »Y2»Z1) 

1929  LET  Y 3  =  Y +  H*H3 

1339  LET  x2=FNG( T  +  H) 

1949  LET  Z2=( FNG( T>3*h/2 ) -FNG( T>H/2 ) ) /H 

1959  LET  K4=FNF(x2»Y3*Z2) 

1  3 62  LET  Y=Y*h*(*W*K2*2*R3  +  R4) /8«<3 
1  2  79  LET  T=  T  +  H 

1990  REM  CHECK  TO  SEE  IF  IT'S  T|m£  TO  PRINT  OUT. 

1999  If  (  T>H-b)*S6N(L)>1  .0E-9>=9  .9  THEN  1120 
1100  LET  H  =  T  H 

1119  GO  TO  1139 
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1  120 

LEI 

R  =o 

1130 

LEI 

A=l NT( E*«)/E+L*(SGN(H>-1.0)/2 

1  1  444 

1  F 

(A+L_h)*$GN(L)“1»3E“9>=0«3  THEN  1160 

1  150 

LET 

A  =A>  L 

1 1  60 

1  F 

( T-A )*SGN( L)> 1 *0E-9>=0 «0  THEN  1190 

1 1  n 

LET 

0  s  A 

1  1d0 

oOSUd  1270 

1190 

1  F 

R  =6  THEN  1220 

1200 

GO 

TO  930 

12  10 

.  REM  AH  E  *E  OGNE  YET? 

1220 

1  F 

Ad5 ( R  -A ) <0  *5E*9  ThEN  1250 

1230 

LET 

0  =R 

1240 

GOSUb  12  70 

12  50 

$  TOP 

12  00 

REM  ITERATION  AND  PRINTOUT  SUBROUTINE 

12  U 

LET  Hi =0-1 

12  00 

LET 

Z 3  =(  FNG<  T  +  H1 /2 ) -PNG(  T-ri 1/2 ) )  /H  1 

1290 

LET 

X  3  =  FNG( T) 

13  0'4 

LET 

R  5  =  F  N  F  ( <3  *f  *Z3) 

13  10 

LET 

Z4s(F\G(T*H1)-F,NG(I))/Hl 

1320 

LET 

x  4  =  FNG( Tf Hi /2 ) 

1330 

LET 

Y4=YhH1*R5/2 

1340 

LET 

K6=pNF(  X4»Y  4»Z4) 

13  50 

LEI 

Y  S=Y+H 1 *R  6/2 

1  3  b0 

LET 

R  7  =  F N F (  X  4  »y  5  >Z4) 

13  70 

LET 

Y  o  =  Y  1 *R  7 

133  0 

LEI 

X  5  =  F  N  G  (  T +  H  1  ) 

1390 

LEI 

ZS*(P.NG<  [♦3*nl/2)-FNG(  T  +  Hl/2)  )/H1 

140  0 

LET 

R  3  =  F  N  r  (  x  5  >  Y  6  *  Z  5  ) 

1410 

LET 

Y  7  =  Y>H1*(R5H2*RbH2*R  7>rt8  )/6.9 

1  420 

1  F 

M=0«0  ThEN  1460 

1430 

LET 

0  =  (  Y  7-Y ) /Hi 

14  40 

K  R  1  N  T  I  *  X  3  >  17*  Z3  *  0 

14  50 

GO 

TO  14  70 

1  4o0 

PRINT  I ♦  X  3  *  Y  7 

1470 

RETURN 

1480 

END 

R  E  A  DY 

81 


L  I  5  I 


bIDUIA  23  :  53  03-DEC-75 


100 ' 

M  '0  6  —  d  1  DL0  1  4 

11k)' 

12/  ' 

UtSCnlPTI  0.N--A  MODIFICATION  OF  H 1 OLtfl  • 

Th  |  S 

h P OGr  A  M 

SOLVES 

13  /  ' 

IHE  INITIAL  VALUE  PROBLEM  DEFINED  FOR  61 

DL0  1  • 

IN  THIS  CAiE  * 

1  4/  ' 

IHE  MAXIMUM  INPUT  VOLTAGE  LEVEL  IS  INCREASED 

from  10 

VOL  TS 

1  50 ' 

10  20  VOLTS  IN  STEPS  OP  0.25  VOLTS.  The 

PEAK 

OUTPUT 

VUL  TA  GE 

1  6/  ' 

IS  PRINTED  ALONG  <*ITH  ThE  CORRESPONDING 

1  NP  U  T 

VOL  TAyE 

• 

1  l<d  ' 

IhtSh  VALUES  PERMIT  A  D  E  f  EP  '*1 1  NATION  OF  LINEARITY  FGk 

IHE 

Id v?  ' 

:>  a  L  E  D  OUTPUI  hEPkESENIED  d  Y  HE  U  U  H  U  f 

h  U  L  S  t  • 

190  ' 

2  d  0  ' 

210' 

220  ' 

230 

K  E  M  INSERT  THE  h  R  1  N  T  0  U  T  TITLE* 

2  40 

PRINT  “PEAK  OUTPUT  DIODE  VOLTAGE  AS  A  FUNCTION 

OF** 

d  54 

Print  ••  peak  channel  voltage" 

2  o0 

h  r\  1  to  T 

2  U 

pRI  NT  “  X  PEAK  "  »“  r  P  E  A  i\  “ 

2.30 

Kc*1  DEFINE  ThE  DlFhtrtEtoMAL  EQUATION* 

290 

OEF  pnp( x  ,r *Z  ) 

30  0 

LET  C=1 .4755E-1 4/S0r( AbS(x-y)) 

3  10 

LET  0  =  2.9 51 E- 15/S  OR ( An S( Y ♦ 2  5-0) ) 

320 

LET  F=C*Z-(2»0E-10)*Y 

330 

LET  G=C-f  0-^  1  •  6 1 8  E  -  1 3 

340 

PNP  =  P/G 

3  50 

P  N  E  N  0 

3  bd 

FOh  1  =10.0  TO  20.0  S  TEP  0.25 

3  70 

REM  INSERT  ITERATION  AND  TABULATION 

DATA. 

38  0 

LET  T0=1.4hb5E-5 

390 

LET  Y  =  -2  4 .0 

400 

LET  d=4«6985E~5 

41  w 

LET  H  =  5 .0E-9 

42  0 

LET  L  =  4. 698  5  E  —  5 

430 

LEI  E=1.0/L 

440 

LET  T=T0 

450 

R EM  DEFINE  THE  DRIVING  FUNCTION 

460 

OEF  FNG(T) 

470 

F  N6  =  -2  5*0 

48  0 

IP  T<15»0E~6  T h 6 N  510 

490 

IF  T>47»0E-b  Th£.'i  510 

00  0 

F  N  0  =  - 1 

510 

FNEND 

520 

IF  (b~T0)*H>0.0  ThEN  540 

530 

LET  H  =  -H 

5  40 

IF  (H-T0)*L>0.0  ThEN  560 

550 

LET  L=-L 

82 


5  o  /  LEI  X  =  F  N 6 ( T ) 

57/  60  TO  7  5/ 

58/  REM  P  E  R  I"  0  rt  M  TnE  6  A  i  |  C  I  T  E  r<  A  T  |  u  Ts  • 

59/  LET  Z=( PNG( T+n/2) -FNG( T-h/2 ) ) /H 

b//  LET  X  =  F NG(  T ) 

bit?  LET  K 1 =FNF( X  »Y  *Z) 

62/  LET  Z1=(  FN6(  T  +  rt) -FNG(  T)  )/n 

63/  LET  Xl=FNG(  T+H/2) 

b4/  LET  Y1=Y+d*Hl/2 

65/  LET  K2=FNF(x1  »Y1  »Z1  ) 

bo/  LET  Y2  =Y  ♦n*r\2 /2 

o  7  /  LET  r\3  =  FNP’(Xl»Y2*Z1) 

od/  LET  Y3=Y«-HM3 

o  9  /  LET  X2  =  F-MG(  N«) 

1  n  A  LET  Z2  =  (Fivj6(T+3*d/2)-ri'«0(T  +  rl/2))/r( 

Jlo  LET  K4sFnf(X2  »Y3  »Z2) 

72/  LET  Y-f*pt*(r\1<-^*r\2  +  2*r'3  +  ^4)/b»^ 

73  o  LET  T=  T  +  ri 

740  G  rt  EC  R  TO  SEE  IF  IT'S  fl^E  T  0  P  R  I  N  T 

750  IF  (  T  +  d-  d)*S0N(L)-*-1./£-9>  =  /’./  T M  E  n  7 a  / 

7  o/  LET  R  =  T  +  H 


77/  GO  TO  79/ 

73/  LET  K=o 

79/  LET  4  =  1  Nl(  E*k)/E*L»(  S  6  N  (  h  )  -1  •  /  )  /  2 

8  /  /  IF  (  A+L-K)*SGN(L)-1  •  /  E  -  9  >  =  /  •  2  T  ri  E  N  82/ 

91/  LET  A  =A*  L 

8/0  IF  (  T-A)*$6N( L)-1 *3S-9>=/ •/  THEN  85/ 
83/  L  E  T  0  =A 

8  4/  GuSUd  9  4/ 

d  5 «/  IF  R  =a  Tn  £  N  3d/ 

8  6/  GO  10  59  / 

97/  REM  ARE  ’*  E  00 N  E  YET? 

dd/  IF  AttS(  H  -A  )  </  »5E-9  Th£i\  91/ 

8  9 v  LET  Q  =  R 

9//  GOSUd  94/ 

9  1/  N  E  a  T  I 


92  Z 
93/ 

9  4/ 

9  5/ 

9  b/ 

9  7/ 
98/ 
99/ 

1  /// 
1/1/ 
1/2/ 


S  TOR 

k £m  ITERATION  AND  PRINTOUT  SUckOuTINE* 
LET  HIsQ-T 

LET  Z3=(FMG( T+Hl/2)-FNG( T-H1/2))/rt1 
LET  X3=FNG(  T) 

LET  K5sFNF(X3»Y»Z3) 

LET  Z4  =  (  FNG(  T*ri1  ) -FNG(  T )  )  /H  1 
LET  X  4=  F  N  G ( T*h1/2) 

LET  Y4=Y+Hl*R5/2 
LET  Kb  =  FNF(X4>Y4»Z4) 

LE  T  Y  5  =  Y  +  H 1 *K  6/2 


OUT. 
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1<Md  LET  ri7=F^F(X4*Yb>Z4) 

1444  LET  Yb=Y+H1**7 
1  v  id  LET  xb  =  Ft\G(T  +  ri1) 

I  y  F  6  LET  Z5  =  (FNG(T  +  3#H1/2)-FNG(T  +  i-I1/2))/H1 
1  7  0  LEI  K8  =  F'\F(Ab*Yb*Z5) 

It)  ay)  LEI  Y7-Y+H1*(rtb+2*ri6+2#K7+r\d)/6*?' 

1494  PRINT  X3  9  Y  7 
1  144  hETUFN 
1114  E  N  0 

*eagy 
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U  3  r 


o  I  OL  d  {  0  1  :  0  2 


0  4- DEC  -  75 


1  00  ' 
1 1  e>  • 
120  ' 
130  ' 

1  40  ' 
1  50  ' 

1  DO  ' 

1  Id  ' 

1.3  0  ' 
Mo  ' 
idd  ' 
210' 
Ltd  ' 

230  ' 

2  40  ' 
2  bw  ' 

2  o0  ' 
2  M  ' 
2b0  ' 


4U  ' 

42  0  ' 

43  1 
4  40  1 
4  D0  1 
4  b/>  1 

4  H  1 
4o  j  1 
4*4  1 
5«cM  1 

5  b)' 
5  ^  ^  1 


NAME--d I DLtf  2 

A  PROGRAM  rO  FIT  A  CURVE  THROUGH  A  S  6  f  OP  DMA 
61  THE  METHOD  Oh  LEAST  SQUARES# 

N  INDICATES  THE  TYpE  OF  CURVt  TO  HE  USED  TO  FIT  TnE  OAfA# 

HE  P  0  *  t  H  FUNCTION  Y  =  A  *  X  T  d  *  I  L  L  rt  E  SELECTED  IF  N  =  <)  • 

HE  E  X  p  0  N  t  N  F I  A  L  FUNCTION  Y=a*EXpH*X)  *  I  L  L  5  E  5  E  L  E  C  T  E  u  IF  \  = 

HE  POLYNOMIAL  Y=0(1)+C(2)»HC(3)**f2*###  +  C(1tf)*XT9  ^  I  L  L  j: 
SELECTED  lh  N  EXCEEDS  1#  (N  WILL  INDICATE  THE  OEGRcE  Or  Hr 

polynomial#  for  example*  n = 3  indicates  a  s  ec  o n  d - o e  or  t e  r u  l  y 

NLHIAL#  NOTE  ThAT  N  IS  HE  N U **o  E r  OF  C  0  n  S  T  A  «\  [  S  IN  Hr 

fglyn04|al--hE\ce*  n  cannot  exceed  ten# 

frit  PrOGrAm  REQUIRES  TrtE  FOLLOWING  I  Nr  U  T  S  • 


290  ' 

570 

LET 

M  =  C  *  **  3 

THE  N U  V rt  E R  Or  DATA  P  A  1  R  d  ♦  NOT  T  u 

30  0  ' 

EXCEED  1 0  0  P  0 1  N  F  d  • 

310' 

5b0 

LE  I 

N  =  £>*3 

TrtE  TYpE  OF  Curve  TO  HE  FITTED  AS 

320  ' 

DISCUSSED  AdOVE#  N  CANNOT  ExCEEu 

330  ' 

1  0  OR  M  WHICHEVER  IS  S  m  A  L  L  t  r  • 

3  40  ' 

59  0 

LET 

1  0=C  *  3 

ThE  DESIRED  OUTPUT  FOFH*  1  Or  2# 

3  50  ' 

60  0 

04  1  4 

X(  1  )  #  Y(1)t  X  H  )  *  Y(2)  Inc  DATA  TO  rt  t  FlTTEu 

3  00  ' 

TO 

1  \  P  4  I  A  3  »  X  4  N  0  r*  4  .3 

3  70  ' 

70  0 

04  T4 

X(M)»  r(«0  MANY  AS  *  ILL  M  T  A  LINE* 

330  ' 

710 

Phil  NT  "THE 

TITLE  OESIPEO  ON  Tic  OuTrUT  P  r,  1  4  i  0  u  [  " 

390  ' 

40  0  ' 

Vi  1  I  rt  1  h 

TnE  P 

R  0  GR  A  M  9 

TnE  FOLLOW  1  NO  0414  lo  C  On  T  4  i  \E  u: 

A  (  I  » J  )  =  TH  E  COEFFICIENTS  OF  ThE  UNKNOWN  CONoTANH 

c ( j )  -TnE  unknown  Constants  that  define  he  fitted  co^ve# 

D(  J )  =  HE  R  I  GH  [  HAND  TERMS  OF  ThE  SIMULTANEOUS  EQUATIONS# 

Ir  10  IS  l#tf*  THE  UMNO*N  CONSTANTS  AND  ThE  DATA  POlNH  AkE 
pnInTED  OUT#  IF  I  IS  2#H  ThE  SYSTEM  OF  LINEAR  EQUATIONS 
IS  PHI  N  T  E  0  *  AS  a  E  L  L  Ad  THE  U  NK  N  0  * (M  CONSTANTS  AND  ThE  DATA 
POINTS#  ThE  S  U  *  Or  THE  SQUARES  Or  THE  ERRORS  IS  A  {_  *  A  Y  d 
p  r  I  n  T  E  D  # 
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ddw  1 

»  *  *  *  *  *  *  *  MAIN  Pr<  QG*  4 

«  *  * 

54  0  • 

d5  c 

Ul  A  A( 1/  >10)  »  C(  10)  *  D( 10) 

*  X  (  1  V  w  )  > 

r  ( it)*/) 

b  o  e? 

K  E  d  rt  EA  U  1  iNiK  U  I  DATA. 

5  In 

LET  d=4 1 

bd  / 

LET  N=4 

59  4 

LEI  10=1 

b/) 

DaTA  It)  *2’  22  *492  8’  1  •  2  5  ^ 

2  2  •  5  r  J  1  vX  9 

1^*5’ 

22 

•  5v>  9  1’  1 1  •  7  b  ’ 

44 

*d1  7j 

bid 

DATA  11.  22  •  52  5r  >  11.25* 

2  2  •  53  3  b  * 

11.5’ 

2  <; 

•5431’  11.7b’ 

22 

•  5  5  /  4 

o  4  *6 

u  -•  T  A  1 2*2*  4  2  •  5  3  d  8  *  1  2  *2  ^  * 

22  •  5f 72  ♦ 

1  /  • :  ’ 

22 

•  5  7  ^  c ’  12*75* 

44 

•  ’0  4* 

8  i  / 

J  a  T  A  13*  2  ’  22*5^25’  13  •  2  5  ^ 

4  2  •  ^  4  1  /*  * 

13*5’ 

22 

•  OOF’  13*  7  c  * 

22 

•  A  1  o  4 

o4  i 

DATA  14*/)*  22  *62  6d’  1  4  •  2  b  » 

22  .-3  55  * 

14.5’ 

22 

•  r.  4  4  3  ’  14*75’ 

22 

•  c  5  3  »' 

b  3*2 

DATA  1  b  •  ^  2  2  •  8 1>  1  d  *  1  5  •  2  5  ^ 

22.6  7/  7  * 

1  5  •  5  * 

22 

•  879b’  15*  7  d ’ 

22 

•  r)doc 

bb/ 

u  a  T  A  16*2  9  22*69  78*  1  b  •  2  5  * 

2  2  •  7  4  b  7  ♦ 

lb. 5’ 

22 

•  7  1  cd  ’  1  D  *  /  D  ’ 

24 

•  72  bt 

0  7/ 

DATA  17*2’  22*7343  *  1  7  •  2  o  ♦ 

2  2  •  7  43  c  ♦ 

17.5’ 

22 

•  7  53  t)  ’  1  7  *  7  b  ’ 

2  4 

•  7  o2  b 

cd  .5 

Data  1d*v/»  22  •  7  72  1  *  13  •  2  b  ? 

2  2  •  7  d  1  3  * 

1 6  •  5  ’ 

22 

♦  7  9  1  o  ’  1  d  •  7  c  * 

44 

•  a  /  1  5 

bd  0 

u  A  T  A  1  d  •  *•  9  22.8115’  19*25’ 

22*3^1c* 

19*5’ 

22 

.1319’  lfi*7c’ 

4  2 

*  d  4c  d 

7  /  d 

*j  A  T  a  2  '4  •  v5  *  2  2  •  B  5  ^  -1 

7  1  0 

r  n  1  0  [  '•  0 U  f  9  U  F  VOLTA  3d  -  > 

A  r*J\C  f  1 

;j  i  vi  u  r  1 

ir 

U  F  VULfAGc'* 

7  A  <■• 

ml  ^  T 

/o/ 

h  un  a  -  1  f  0  * 

7  c, 

r<  EA  U  X  (  A  ) 

7  dt 

Kcu  r  (  a  ) 

7  c/ 

id  £  A  I  A 

7  7/ 

h  uAlCULA  i£  L  0  *j  A  t<  1  I  r)  i*i  ^  J  I-  1  vi  H  U  f  m  D  1  < 

AT£  DATA  |i- 

7o  i 

r\  £  >i  ^  .>  £  a  La  £  X  r  U  f\  £  '  I 

1  A  L  C  Urt  V  E 

1  o  -El 

•*  vC 

H  T  T  E  u  • 

7  j  .j 

1  h  *\i  >  =  2  •  /  I  h  £  ;  37  <:• 

o  «5  / 

h  jn  ^  =  1  F  0  *• 

d  1C 

1  r  T  (  a  )  >  /  •  /  M  o  5  *i 

12/ 

KKlNf  ,,\E‘jAllVc  VALucd  r 

C«  r:  C  A 

Ni\Ol  Ob 

c 

A  r  U  A  P  A  Nj  .  ^  1  1 

<  • 

** 

S  w  ^ 

r'  K  1  N  i  *4  A  •  i  £  X  r  0  N  £  IN  F  !  A  L  r  1 

[ 

0  4  /  5  f  Ur 

o  5  c'  LEI  Y(ft)=LUii(f(^I) 

dot'  <  t  <  I  b 

b  7  /  r<  E  X  CiLLuL-ilE  L0>jAr!lfr-!-'^  Of  I'if'Uf  i^J1,  i  so  j  l  c  I  i 

d  d  /  c  •’  Id  4  r  v  >•  t~  u  jr  V  c  I  i  o  f  I  i  o  r  I  I  I  c  U  • 

o  j  ■/>  If  u  -  1  •  •>’  I  m  E  in  d  7  / 

d  .  riir!  r\=1  I • 
d  1  2  lr  *  (  b  )  >  £  •  /  I  d  r  ,  j  4,1 

d2«  r  -  |  (r  •*  i\iEG4  T  I  VE  VALUES  FOR  x:  CANNOT  USE  A  ro  f  !  I  u  -  r  I  f  - 

d  3  d  5  T  OF 

d'i  >  LEI  aU)=L0G(X(K)  ) 

i 50  NEXT  K 

)  r>  /  RE'*  UALCSJLAIE  I  r|  t  EL  t'lcN  I  5  Of  I  n  i  A  -  A  I  Pt  I  /  A  •;  L  I  rl  C  u*/:^  lu'  • 

d  7 0  LEI  N 1 =  N 

d  o  x  If  \  1  >  -  2  »  '  r  "t  E  N  1  iNix 

did  LEI  N1=2«0 

1  <■  .)  r  Urt  1=1  I  0  N  1 

1  /•  1  &  rOrt  J  =  1  10  N1 

\Uc  If  (  I  +J)  >2  ••?  Then  tv  5vd 

1  /  3  *  L  E  I  A  (  I  ♦  J  )  =  v, 

10  4/  GO  10  1/1  dv! 

1/5.)  LEI  A  (  I  *  J  )  =  /  « •? 
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Mb/  rOR  *  =  1  10  'I 

\rl<)  LET  6(  |  #J)=A(  I  »J)  +  *(r\)t(  MJ-2) 

1  ft  8  /  IN  E  A  I  R 
I  /  i  /?  IM  P  A  [  J 

I  1*M  LET  D(  I  )=/•/ 

1111/)  FOR  K  -  1  TO  M 

I I  2C?  I  F  I  >1.M  i i E \  115/ 

1 1  3 w  LET  D(  t  )  =D(  IMY(a) 

114/  GG  TO  116/ 

115/  LET  D(  I  ) =0(  I  )+Y(*)*<(4)t(  I  -1  ) 

1  1  n  /  |\  E  x  r  K 
1  1  7/  IN  Ex T  I 

113/  REm  WRITE  ME  Oh  T I  OMAL  OUTPUT  IF  M=2* 

113/  IF  I  /  <  2  •  /  T  ri  E  ^  13// 

12  0/  PrI  NT  “  C  0  E  F  F  I  C  I  ENTS  IN  ME  a  Y  $  T  E  *  OF  LINEAR  EO  UA  T  I  ur.s  :  - 

I  ^  1  ✓.«  F  r<  I  NT 

122/  M A  T  PRINT  4 
1  2 .3  v  PRINT 

12  4/  K  ft  I  N  T  '*  M  r  ft  I  G  ri  T  n  A  .  i  u  C  0  L  0  -  ’  <  VECTOR  IS  .* 

1  2  S ft  PRINT 
1 2  o/  FOR  I  =1  TO  N1 
1  2  7  v)  P  R  I  N  T  D  (  I  ) 

12  b/  h  E  x T  I 

12//)  rt  EM  SOLVE  M£  SYSTEM  OF  SIMULTANEOUS  LINEAR  EQUATIONS* 

ISO/  GO  o  U o  13  y 0 

1314  K  EM  wRIIE  Me  OUTPUT* 

13  2/  F  ft  I  IN  T 

133/  I  F  N >  =  1  ♦  ')  MEN  13  hi 
1340  LEI  C  v)  =  E  A  p  (  C  (  1)  ) 

13  5/  PRINT  ••pOwEp  FUNCTION  CUftvt  FIT:  Y  =  ■*  ;  Cft ;  *  a  T“  ;  C  (  2  ) 

1  3  btf  GO  TO  1  53  4 

13  70  IF  iN >  1  •  /  MEN  141/ 

138/  LET  C  /  =  E  <  R  (  C  (  1  )  ) 

134/  r  r  I  NT  41  £a  ft  u  \  E  N  T  I  A  L  C  Jr  v  £  FIT:  T=“;C^i"*E^t"iC(2)U'*0" 

1  4^i /  o J  TO  153/ 

14  14  I  r  \i  =  1  /  •  /  1  fi  E  N  14  E  A 
142/  LET  n2  =  M  1 

14j,;  r  0  R  I  -  *  c  \  0  I  ft 
14  4/  LET  C  (  I  )=/’•/. 

1  4  5v)  "j  E  ^  T  I 


1  4  bid 

PHI  i'l  r 

"  p u  L  Y  m  0  m 1 AL 

C'J-i  Vc  PIT: 

" 

WO 

Pn  1  NT 

1  Jrsir'* 

rrJ  \T 

"  r  s  (  •' :  C  (  1  )  :  '• 

)  *  (  ;  0  ■  4  )  : 

"  )  •  \  «•  ( •• ;  C  (  j ) ; 

1  49  " 

P  rt  1  i'll  I 

il 

•*■(">0(4); 

•• )  *<  T3  +  (  "  ;0(  d)  : 

14  )  *  A  T  d+  •  •  •  *' 

1 

ph  1  N  f 

II 

■*•(":  C  (  R ) ; 

'• ;  *  x  r  5-r  ( ;  0  (  / ) ; 

44  )  *  *  f  •  *  •  •* 

151-* 

PHI  NT 

i« 

*  ( " ;  C  (  h  ) ; 

•,)*<r 7  +  ( •• ;  ;(•>); 

44  )  *X  Tb  +  •  •  •“ 

1  52  it 

PH  1  NT 

ll 

*  ( " ;  c  ( i  i-i ) 

;  '• ) *x  r°" 

1  53  y 

Prt  1  N  T 

1  54«; 

H  Erf 

*hIT£  OUT  Tit  InpiJT  VALUES  Or  <  AMO 

M  A  N  C  T  h  E 

1  5  b  nJ 

REM 

C  OK  RES  PON  01 

NO  CALCULI  TED  VALUES  Or  i 

• 

1  5  oU 

1  F  N  > 

=2  Th£N  1  o 4 » 

1  5  /  Id 

pOh  *  =  t 

TO  M 

158/  LET  r(r\)=  E^P(YCK)  } 
1  59 NEXT  K 
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1  0  3  /J 

1614) 
1  fil'd 
1  b3  3 
1643 
1  b50 
1  bb3 
1  b73 
1  66  3 
168'/) 
1  73  2 
171/» 
1  722 
1  73  * 
1  74/1 
1  7  5  2 

17  6k) 
1  7  73 
1  7b  2 
1  79'/5 
1  b  /}  3 
1b  l/i 
1  b  a  3 
1  b  3  v) 

18  4k) 
1  0  D/‘ 
1 8  o/) 
1o  73 
I860 
1  bS/, 
19/)  2 
1914J 
1934 
193/) 
194/) 


IF  N  =  1  •  2  THEN  1  b  4  ^ 
r iJK  6  =  1  TO  "i 

LET  X(6)=EXR( <(6) ) 

.9 Ex :  <\ 

bRINT  '‘Dili.  POINTS  AND  JALOJLATEJ  VALUES  Ob  Y" 

P  tk  I  1/  f 

P  R  I  f4  r  “VALUE  OF-  X“>-‘ VALUE  OF  r  “  *  "  C  A  LC  UL  A  I  E  D  Y  " 
PR  I  r«  T 

LET  33  =/<  •  £ 
rOri  6  =  1  T'J  9 

I  F  N>  =  1  .3  Hen  1  73 1 
LEI  Z  4  =  C  2*  v  ( a  )  tC  (  2  ) 

Ou  fu  1b U 

IF  \  >  1  •  3  7  ri  E  ■>i  17b  f; 

LET  Z 4  =  0 2* E‘ P ( C ( 2 )  •  X ( a  )  ) 

GO  10  16  1*' 

LET  Z=C(  1  )*C(2)*X(a)+C(3)  *<(  a)  12 
LEI  Z1=C(4)*/(b)  *3*CLS)*x(r')  f4 
LET  Z2  =  0(o)*<(rt)l5+0(7)*A('')rb 
LET  Z  3  =  C  (  6  )  *  X  (  a  )  t7-»C(o)»x(k)  To 
LET  Z  Y  =  Z>Z  1 +  Z  2  v  Z  3  +  C  (  U')*x(r.)  T9 
LET  Sv:  =  cj;/*(Z4-t(i\))t2 
r'H  I N  T  X  (  A  )  *  Y  (  3  )  t  Z  l 
w6X  T  a 

R  E  M  n  R  I  T  E  T  b  E  E  a  k  C  R  6  E  T  «  E  c  \  I  N  a  U  T  A  1 0  Z  A  L  0  U  l 
Pp.  I  lil 

PRINT  “THE  5UW  OF  T  i-t  £  SO  uAaE  OF  TbE  Ebb  OR  3  =  " ; 

3  T OP 

REM  SUBROUTINE  TO  SOLVE  SIMULTANEOUS  LINEAR 
FOR  1=1  TO  N 1 
FOH  K  =  1  TO  N 1 

IF  6  =  1  f  h  E  N  2 /dd 
LET  C3  =  -A(K  >  I )  / A (  I  »  I  ) 

FOR  J  =  1  TO  N  1 

LET  A(a/J)=A(6  »J)*C/i*A(  I  >  J ) 


1  9  Ski  lb  J  =  1  *  vX  TnEN  19  73 

1 9 b4i  GO  TO  1  *c  3 

1973  LET  A(K>J)=3.3 
I960  NEX T  J 

199  0  LET  i)(  h  )  =D(  A  )  >C3*  U(  I  ) 

2  3  t;  3  h  E  X  T  R 

2313  LET  Co  =  A ( I  / 1  ) 

2323  FOR  J=1  TO  91 

2333  LET  A ( I  » J ) =  A (  I  •J)/C3 

2  3 43  NEXT  J 

2  3  62  LET  A  (  I  » I  )  =  1  •  3 

2  3  b3  LET  0(  I  )=0(  I  )  /C3 

2  3  7w  N  £  x  T  I 


2363  PUR  1=1  TO  N 1 
2393  LEI  C (  I  )  =D(  I  ) 
2  133.  NEX  T  I 
2113  RETURN 
2  12  3  END 


A  T  E !.  b  0 1  N  i  S  • 
S3 

EO  UA  I  I  0  ii  « 
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L  I  5  T 


5  I  JL-4  j  2  3  •  1  4 


,•‘4-060 
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1  on  '  '■» E  -  -d  I  DL/  3 

1 1  w  ' 

12/)'  CESCH I  P  T I  ON  — A  PkOGbA.i  TO  EVALUATE  THE  G'jTpUT  VOLTAGE 

13/'  oh  A  0  S  F  E  T  SOURCE  FCLLO'/iE?  /» 1  T  H  T  n  E  INPUT  ppOm  A 

14.!'  CCO  CHANNEL  THROUGH  a  hEVEhSE  tIASEC  DIODE.  T  a  fc  CCD 

15/)'  channel  voltage  mi  i  e  x  ♦  the  input  to  the  source 

1b/1'  i-OLLO«E~.  a«  T  L  L  hE  v  1  »  HE  0  U  T  h  J  T  Oh  HE  S  f 1 U  h  0  E  r  C  L  L  0  a  E  h 

17/)'  IS  r  . 

lb,)  ' 

1-o'  THE  V.05FET  I  Np  !J  T  VOLTAGE  IS  Appw  Ox  I  -<tA  TED  -  r  Hr 
2/)2'  POLY. NO -I  AL  RELATIONSHIP: 

2  1  <■'  ' 

22X'  Vl  =  -(A  H-W  +  CHX  f2-'j1«*  T  ?  ) 

2  3/  ' 


2  4  y)  1 

1  HE 

.•*  0  S  r  £  I  DR  A  1  \  -  o  'j  ;J-h  C  E  C  J,vr(  ;  \  I  IS  A  r 

ph Ox  I  .  a  TEC  h  r 

2  52  1 

1 

2  5  /*  ' 

'  1  = 

-/*(  VI -Y -V/  )  f 2 

2  7/5  1 

1 

2  i  /  1 

1  *hP  n 

P: 

2  2 

1 

r\  -  THr  GAIN  r  A  .«<  A  V  P  T  E  IN  A  E*  £  /  VG  L  f 

5  /)  v) 

1 

V  ^  r  TmE  THkES^uLu  voltage  in 

V  l  i  L  T  5 

31 1  1 

1 

3  2/5  ' 

1  THE 

SOURCE  KOI LOrfER  LOAD  RESISTANCE  1 

S  R  • 

3  3  1 

1 

0  4/  ' 

3  52  1 

1  *  * 

3  6/  1 

1 

4  70 

P  N  1  N  T 

"OUTPUT  Oh  SOUrtCE  FOLLU.vEh  AS  A 

E  U  hC  f  1  U  N  C  r 11 

3s/) 

pp  1  NT 

••  Channel  voltage  in  CCD" 

34/ 

pH  1  NT 

43  0 

P  H  1  N  T 

"CHANNEL"  » "  0  U  T  P IJ  T  IN" 

41/3 

HK  1  NT 

" VOL  TAGS" »" VOL TAGE" »"SA  PJP A  Tl UN" 

42/ 

PHI  NT 

43  3 

LET 

K=2  .2/32E-6 

4  4/ 

LET 

V/ =-  1  .  76465 

4  5/ 

LET 

P  =  5  •  /  E  4 

4  5/ 

LET 

41=22.15/ 

4  7/ 

LET 

r)1=3  *8/4  4 1 -2 

48/ 

LET 

CT  =  -n.i51')5E-4 

44/ 

LET 

01=2  .33133E-5 

5  /  3 

LET 

• 

J"> 

'■VI 

1 

If 

UJ 

31/ 

F  OH  X 

=1/  10  2/  STEp  75  •  2 5 

52  w 

LET 

C2=01*x+C 1 

53/ 

LET 

b2  =C2*x  +  h 1 

54/ 

LET 

VI =h2*  x*A 1 

55/ 

LET 

< 
_ i 

II 

1 

<T 

89 


5b*J 

LEI 

i  =  2  •  (  v  a  -  \n  )  + 

5  h) 

LET 

C  =  vH  T2  +  Vt'  t2  - 

5d 

LEI  X 1 = - x 

59 

LET 

r  =SG«( d t2 -4* 

b^b 

LEI 

r  =  (  -3+y  ) /2  •  *> 

0  1  v' 

LEI 

L  =  A S  s (  VI  - r -  V 

Od/* 

LEI 

J=AoS(E-Y)  - 

63  ^ 

1  F  L< 

= J  ThEN  6  c  d 

64^ 

K  F  1  i'i  T 

a  1  »  y  *  "  N  0  “ 

6  5k> 

GO  TO 

67d 

6b^ 

r  v,  1  N  T 

xi*  r  *  "yes 

b  70 

NtX  T  X 

bdtf 

S  TOF 

69*J 

fc  NO 

n  EA  uY 
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